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Abstract 

Josephson-j unction arrays are ideal model systems to study a variety of phe- 
nomena such as phase transitions, frustration effects, vortex dynamics and 
chaos. In this review, we focus on the quantum dynamical properties of low- 
capacitance Josephson-j unction arrays. The two characteristic energy scales 
in these systems are the Josephson energy, associated with the tunneling of 
Cooper pairs between neighboring islands, and the charging energy, which 
is the energy needed to add an extra electron charge to a neutral island. 
The phenomena described in this review stem from the competition between 
single-electron effects with the Josephson effect. They give rise to (quantum) 
Superconductor-Insulator phase transitions that occur when the ratio between 
the coupling constants is varied or when the external fields are varied. We de- 
scribe the dependence of the various control parameters on the phase diagram 
and the transport properties close to the quantum critical points. On the su- 
perconducting side of the transition, vortices are the topological excitations. 
In low-capacitance junction arrays these vortices behave as massive parti- 
cles that exhibit quantum behavior. We review the various quantum-vortex 
experiments and theoretical treatments of their quantum dynamics. 
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I. INTRODUCTION 



A. Josephson-junction arrays 

The first artificially fabricated Josephson- Junctions Arrays (JJAs) were realized twenty 
years ago at IBM [|T| as part of their effort to develop an electronics based on superconducting 
devices. In the first ten years of their existence, Josephson arrays were intensively studied 
to explore a wealth of classical phenomena . JJAs proved to be an ideal model system 
in which classical phase transitions, frustration effects, classical vortex dynamics, non-linear 
dynamics and chaos could be studied in a controlled way. The observation of the Berezinskii- 
Kosterlitz-Thouless (BKT) transition in Josephson arrays is probably one of the most 
spectacular experiments |§] in this respect. 

All classical phenomena can be successfully explained by studying the classical 
(thermo) dynamics of the phases of the superconducting order parameter on each island. 
This approach is justified because experiments are usually carried out at temperatures well 
below the BCS transition temperature. Each island is then superconducting with a well de- 
fined gap, but phase fluctuations are still allowed. Under these conditions, classical JJAs are 
a physical realization of the two-dimensional XY-model and above the BKT transition tem- 
perature, phase fluctuations destroy global phase coherence preventing the system to reach 
the superconducting state. Global phase coherence is only restored below temperatures cor- 
responding to the Josephson coupling energy Ej, which is the energy scale associated with 
Cooper pair tunneling between neighboring islands. 

As we now understand, the relatively large junctions at that time had resistances too 
low [| to observe clear quantum effects. By the end of the eighties semiconductor technology 
had pushed device dimensions well below the micron size. It became possible to fabricate 
arrays with Josephson tunnel junctions of sizes 100 x 100 nm^. Circuits with such small 
junctions showed single-electron effects when cooled down to temperatures corresponding 
to the charging energy Ec, the energy needed to add an extra electron charge to a neutral 
island. It was soon realized that the competition between single-electron effects [0,0] and 
the Josephson effect would lead to new, exciting physics. 

An appealing feature of JJAs already emerges at this stage as they can be visualized 
as model systems to investigate quantum (zero-temperature) phase transitions |jl2|,|13[. In 
recent years, this field of research has attracted the attention of many physicists. Experi- 
ments on thin, superconducting films, high-temperature superconductors, spin systems and 
two-dimensional electron gases have all shown the existence of quantum critical points. In 
arrays made of submicron junctions, the quantum fluctuations drive the system through a 
variety of quantum phase transitions. A quantum JJA may be insulating at zero tempera- 
ture even though each island is still superconducting. In the classical limit Ej ^ Ec, the 
system turns superconducting at low temperatures since the fluctuations of the phases are 



^ A simple estimate for value of the junction resistance above which clear quantum effects become 
visible can be obtained by using the Heisenberg relation AEAt > h. By taking the charging energy 
(e^/2C) for AE and r of the order of the junction i?C-time one finds that the junction resistance 
Rn should satisfy the inequality Rn > Rq = /i/4e^ for quantum effects to be observable. 
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weak and the system is globally phase coherent. In the opposite limit, Ej <^ Ec, the array 
becomes a Mott insulator since the charges on each islands are localized and an activation 
energy of the order of Ec is required to transport charges through the system (Coulomb 
blockade of Cooper pairs). Strong quantum fluctuations of the phases prevent the system 
from reaching long-range phase coherence in this regime. 

Granular superconducting thin films are closely related to arrays. In granular films, 
superconducting islands of various sizes and with various coupling energies are connected 
together. Therefore, disorder plays a crucial role in these granular materials while it is 
virtually absent in JJAs (or it can be introduced in a controlled way). Models based on the 
behavior of Josephson arrays also form the starting point to describe the physics of ultra- 
thin, amorphous superconducting films in which superconductivity is quenched by disorder 
or by an applied magnetic field. In these two-dimensional homogeneous films it is believed 
that, although the order parameter is suppressed, phase fluctuations are still responsible for 
driving the system through the Superconductor- Insulator (S-I) transition [Q. 

Another important field of investigation addressed with JJAs, is the study of the quantum 
dynamics of macroscopic objects. In the classical limit vortices are the topological excitations 
that determine the (thermo) dynamic properties of JJAs. In the opposite situation [Ej <^ 
Ec) the charges on each island are the relevant degrees of freedom. Vortices and charges 
play a dual role and many features of JJAs can be observed in the two limits if the role 
of charges and vortices are interchanged. By fabricating arrays with different geometries, 
vortices can be manipulated to a great extent. Quantum dynamics of macroscopic objects 
requires knowledge of the coupling to the surrounding environment [|l^]. To a certain degree, 
the dissipative environment can be modeled and therefore JJAs are prototype systems to 
study macroscopic quantum mechanics as well. Born as a problem related to the foundations 
of quantum mechanics, macroscopic coherence in superconducting nanocircuits is acquiring 
increasing attention since the advent of quantum computation. 



B. Phase-number relation 

Throughout this review, the interplay between the phase of an island and number 
of charge carriers Q on it plays a crucial role. Together they determine the properties of 
quantum Josephson networks. The competition between these two canonically conjugated 
variables is captured by the following Heisenberg relation [p^ : 

[4>hQj] = 2ei 6ij 

where the subscripts i and j label the island positions. 

An elegant illustration of this competition is presented by what became known as the 



Heisenberg transistor |T7| , p!8| . The aim of the experiment was to control and measure the 
quantum phase fluctuations through a modulation of the critical current of the system. 
In Fig.|l| the layout of the device is shown. Two junctions in series (indicated by crosses) 
are connected to a current source. The junction parameters are such that Ej ~ Ec, i.e., 
quantum mechanical fluctuations of the number of Cooper pairs and of the phase of the 
central island are comparable. A large superconducting reservoir is coupled to the island 
through a Superconducting QUantum Interference Device (SQUID). 
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In the experiment the critical current was measured as a function of the apphed flux 
through the SQUID ring. It shows a periodic modulation with a period equal to the su- 
perconducting flux quantum ($o = hjle). The role of the SQUID is to provide a tunable 
coupling to the reservoir of Cooper pairs. When the flux equals an integer times half a flux 
quantum {nPhiQ/2) the coupling is turned off and fluctuations in the number of Cooper 
pairs are suppressed. At the same time, phase fluctuations reach their maximum as indi- 
cated by the Heisenberg relation. At fields equal to zero or an integer number of flux quanta, 
the coupling is maximum so that the amount of charge fluctuations reaches a maximum as 
well. In the experiment, the critical current probes the scale of charge fluctuations: The 
situation with large charge fluctuations corresponds to favorable Cooper-pair tunneling and 
a high critical current. Thus, for zero applied field a high critical current is measured be- 
cause charge fluctuations are at their maximum. At half a flux quantum applied, the critical 
current reaches its minimum value because phase fluctuations are at their maximum. 



C. Structure of the review 

This review is organized as follows. In Chapter the basic physical properties and models 
are introduced. Some theoretical tools to study the phase diagram including the boundary 
of the S-I transition are briefly discussed: the mean-field approximation, the coarse-graining 
approach to derive a Ginzburg-Landau effective free energy, and the Villain transformation 
that leads to a description in terms of charges and vortices. These approaches capture most 



of the essential physics. Sections II C and II D are devoted to a description of the phase 



diagram when charge and/or magnetic frustration are included. Since the number of control 
parameters that can be varied is large, the phase diagram is discussed in some limiting cases 
only. Section [II h| describes the various sources of dissipation in JJAs and their effect on 
the phase diagram. The final three sections of Chapter H report on the transport properties 
close to the S-I transition, the S-I transition in one- dimensional arrays and the physics of 
the field-tuned S-I transitions. In all sections, comparison of the theoretical phase diagrams 
with the experimental results is discussed as well. 

Chapter ^ deals with quantum vortex dynamics. After introducing the important vortex 
properties (vortex mass, pinning potential,...) and its classical equation of motion, a theoret- 
ical description of quantum corrections to the classical equation of motion is presented. The 
remainder of the chapter concerns the description of quantum vortex-experiments. We start 
with the single vortex properties (tunneling, interference, and Bloch oscillations) and then 
proceed with collective vortex motion in quasi-one dimensional samples (Mott insulation of 
vortices and Anderson vortex localization). 

In the last chapter, some future directions are explored. The theoretical description of 
two new experiments is outlined: persistent vortex currents in Corbino geometries and the 
Quantum Hall effect for vortices/charges. The experimental technicalities for the observation 
of these phenomena are described as well. We end this review with a brief discussion 
on Josephson qubits in which fundamental aspects of quantum mechanics and quantum 
information theory can be studied. 

We tried to keep this review self-contained and, at the same time, to give a comprehensive 
overview of the quantum properties of Josephson networks. For each of the sections, we 
present the main ideas without going into a detailed enumeration of all the results obtained in 
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the field. There are some topics which are not discussed here. Probably the most important, 
which would require a review by itself, is the effect of disorder which seems to be more 
important for granular materials and ultra-thin films. 

Basics in superconductivity and Josephson physics can be found in the books by Tin- 
kham [|1^] and by Barone and Paterno [^. Since many ideas (e.g. persistent current, 
localization) were born in the field of more traditional mesoscopic physics we refer for these 
topics to the books by Beenakker PT| and Imry [E^] and to the conference proceedings PB|,E3 



Various other aspects of JJAs have already been discussed in previous reviews devoted to 
this topic 

Throughout the review we put h = ks = c = 1. Distances are expressed in units of the 
lattice constant a. We restore S.I. units in the formulas expressing measurable quantities. 



II. QUANTUM PHASE TRANSITIONS 

A quantum Josephson array consists of a regular network of superconducting islands 
weakly coupled by tunnel junctions. Thanks to submicron lithography, array's parameters 
(associated to the shape of the islands, the thickness of the oxide barrier,...) can be made 
uniform (virtually identical) across the whole array. With present-day technology variations 
in junction parameters are below 20 % across the array. The dimensions of the unit cell are 
of the other of few /im^ while the superconducting islands have an area of about 1 fim^. The 
largest samples consist of about 10000 junctions (e.g. 100 by 100 or 1000 by 7). 

Quantum arrays are fabricated of all-aluminum high-quality Josephson tunnel junctions 
with a shadow-evaporation technique. The evaporation mask is made of electron sensitive 
resists in which the bottom resist layer has an undercut to ensure a proper lift-off after 
evaporation. Junctions are formed by evaporating a thin aluminum layer (25 nm) as the 
bottom electrode followed by in situ oxygen oxidation and evaporation of the counter elec- 
trode of about 50 nm from an opposite angle. A more detailed discussion of the fabrication 
techniques is presented in the Appendix ^ 



A. The model of a Josephson Junction Array 

1. Quantum Phase Model 

In Fig. 1^, we show a scanning-electron-microscope (SEM) picture of a Josephson junction 
array. Its schematic representation is given in Fig. ^. In this square geometry the coordina- 
tion number, z = 4. The coupling strength between adjacent islands is determined by the 
Josephson energy Ej = $o-^c/(27r). This coupling energy is inversely proportional to the 
normal-state junction resistance R^. Experimentally Rn is determined from the normal- 
state array resistance measured at 4.2 K, R]\r = {N + l)r]yMy/Mx, where My is the 
number of cells across the array and is the number of cells along its length. The maxi- 
mum junction critical current Jc, in the absence of charging effects and thermal fluctuations, 
is assumed to be given by the Ambegaokar-Baratoff value |^ 

A 

IcRn = tt— 
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with the measured critical temperature Tc- For a BCS critical temperature Tc = 1.35 K one 
gets Jc-RAr=0.32 mV at low temperatures. 

Quantum effects in Josephson arrays come into play when the charging energy (associated 
with non-neutral charge configurations of the islands) is comparable with the Josephson 
coupling (the physics associated with charging effects in single normal and superconducting 



junctions has been reviewed in Refs. [|TO|,|Tll). In addition, as explained in the introduction, 



the junction resistance should be of the order of (or larger than) Rq |^ . 

Arrays are made in a planar geometry, in which each island is coupled to each of the other 
islands and to a far-away ground by its self-capacitance Cq. The junctions are made of two 
overlapping superconducting layers separated by a thin oxide layer and the main contribution 
to the capacitance therefore comes from the junction capacitance C. An estimate of the total 
island capacitance Cs is obtained from measuring the voltage offset (K//set) in the I - V 
characteristics at high bias currents at T=10mK in a magnetic field of 2T. Using the so- 
called local rule [PI1| , PD| ] , Cy, = Me^ /2Voffset- For junctions of 0.01 /im^, Cy, is found to be 
1.1 fF. If one identifies C with C^, the specific capacitance is 110 fF //im^. Measurements 
on large- area junctions have yielded a specific capacitance that is about a factor of two 
lower. This discrepancy shows that stray capacitance (capacitance between next-nearest 
and further neighbors) may play a role in Josephson circuits as pointed out by Lu et al. pi 



However, for simplicity one often identifies the measured as the junction capacitance C. 

Reliable estimates of this self-capacitance (Co)are obtained from separate measurements 
on small series arrays with high Ec/Ej ratio. A magnetic field of 2 Tesla is applied so that 
the series arrays are in the normal state; Co is then measured by varying the potential of 
the circuit with respect to the ground potential. Recording the current through the circuit 
yields a periodic signal with period c/Cq. For islands of 1 /im by 1 /xm, Cq ~ 1.2x10"^^ F 
which is much smaller than C 

As already mentioned, the electrostatic energy can be determined once the capacitance 
matrix Cij and the gate voltages (if present) are known [^,^. Generally one only con- 



siders the junction capacitance C and the capacitance to the ground Cq. In this case the 
capacitance matrix has the form Ca = Cq + zC, Cij = —C (if i,j nearest neighbors) and 
zero in all other cases. Consequently the charging energy (for two charges placed in islands 
i and j of coordinates and respectively) is given by 



,2 .2 . gik.(r,-r,) 

2^'^ 2 J 47r2 Co + 2C(1 - cos k^) +2C{1- cos ky) 



which is well approximated by the expression: 

Elf . ^K„ (^) (2) 

Here, Kq is the modified Bessel function. The charging interaction increases logarithmically 
up to distances of the order of the screening length A and then dies out exponentially. The 
characteristic energy scale is 

Er< = 

^ 2C 
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Equation (0) assumes three-dimensional screening and the range of the electrostatic inter- 
action between Cooper pairs is given by ( in units of the lattice spacing): 

X = ^/c/c'o . 

. If the two-dimensional limit is considered (if e.g. the array is sandwiched between two me- 
dia with large dielectric coefficients) the screening length scales linearly with C/Cq yielding 
a longer ranged interaction. 

From all these considerations, we arrive at the following Hamiltonian which describes 
Cooper pair tunneling in superconducting quantum networks (quasi-particle tunneling is 
ignored at this stage). This model is frequently called the Quantum Phase Model (QPM) 
and is its most general form it is given by: 

H = Hch + Hj 

= I Y.iQ^ ~ Qx,,) Cr^ {Q, - Q.,j) -EjY. cos (0,, - 0, - A,,) . (3) 

i,j <i,j> 

The first term in the Hamiltonian is the charging energy in which the C^j^ is the capacitance 
matrix; the second is due to the Josephson tunneling. An external gate voltage K,i gives the 
contribution to the energy via the induced charge Qx,i = 2egx = J2jCijVx,y This external 
voltage can be either applied to the ground plane or it may be (unintentionally) induced by 
trapped charges in the substrate. In this latter case Qx,i will be a random variable. 

A perpendicular magnetic field with vector potential A enters the Hamiltonian of Eq. (^ 
through Aij = 2e // A-dl. The relevant parameter that describes the magnetic frustration is 

/=(l/27r) 

p 

where the summation runs over an elementary plaquette. In quantum arrays, the 2D fiux 
penetration depth A^(T) = $o/2vr/io/c(T) is much larger than the array size so that the 
magnetic field is essentially uniform over the whole array, i.e., / is position independent. A 
similar conclusion can be drawn by considering the ratio of the geometric inductance (we 
estimate it to be of the order of 1 pH) to the Josephson inductance (larger than 1 nH). 

Throughout this review, two limits of the QPM are frequently discussed: C 3> Co and 
C <^ Cq. The former limit has already been discussed in detail as it is the appropriate 
regime of Josephson arrays. The latter limit is more appropriate for granular films that 
have a short-range Coulomb interaction. To describe these systems we use the following 
notation. When the on-site contribution is dominant, the characteristic energy is 

e2 _^ 

However, some properties (see e.g. discussion in Section PD| ) are crucially dependent on 
the details of the electrostatic energy at small distances, i.e., on weather the nearest-neigbor 
interaction is also included or not. Ei represents the variable denoting this nearest-neigbor 
interaction; E2 the interaction between next nearest-neighors and so on. 

The two contributions in the Hamiltonian of Eq.(^) favor different types of ground states. 
The Josephson energy tends to establish phase coherence which can be achieved if super- 
currents fiow through the array. On the other hand the charging energy favors charge 
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localization on each island and therefore tends to suppress superconducting coherence. This 
interplay becomes evident if one recalls the Josephson relation (which here can be obtained 
at the operator level by calculating the Heisenberg equation of motion for the phase) 



_ 2e _ 2e 1 
dt h ' h 



= i-cr,'Qj (4) 



A constant (in time) charge on the islands implies strong fluctuations in the phases. On the 
other hand phase coherence leads to strong fluctuations in the charge. 

The low-lying excitations of the model defined in Eq.(^), are long wave-length phase 
waves whose dispersion relation can be obtained by considering the QPM in the harmonic 
approximation 



i,3 <i,j> 



The dispersion relation of these modes (usually named spin- waves from the magnetic analogy 
of the Josephson coupling with the XY model) depends on the form of the capacitance matrix 
(see Section pIB| ). The QPM possesses topological excitations as well, charges and vortices, 



that will be discussed in Section II B 3 



A qualitative understanding of the phase diagram can be obtained by considering the 
two limiting cases in which one of the two coupling energies is largest. For simplicity we 
look at the ground state of the system ignoring external voltages and magnetic field. If the 
Josephson term is dominant, the array minimizes its energy by aligning all the phases, i.e. 
it is in the superconducting state. If instead the charging energy is dominant, each island 
has a zero charge in the ground state. In order to put an extra charge on the island one has 
to overcome a Coulomb gap of the order of the charging energy (max Ec,Eq). The array 
behaves as an insulator although each island is still in the superconducting state. 



2. Dissipative Models 



Since the seminal paper by Caldeira and Leggett p6| , |37[] it became clear that dissipation 
changes the quantum dynamics of macroscopic systems. One can therefore ask the question 
to what extent dissipation plays a role in Josepson-junction arrays and what its role is 
on the SI transition in quantum arrays. At low temperatures one expects quasi-particle 
tunneling not to be present since the charging energy is smaller than the superconducting 
gap A. Experiments on small arrays indicate that even at mK-temperatures a small but 
finite amount of quasi-particles is always present, although it has not been possible to 
discriminate the exact details of dissipation. Therefore, we treat the various models that 
have been proposed to describe dissipation in superconducting networks. 

The QPM defined in Eq. (^) only accounts for Cooper pair tunneling between neighboring 
islands and needs to be generalized. The appropriate description is formulated in terms of 
an effective action by the authors of Ref . fT^JS^ : 

Z = JllD<l>,{T)exp[-S{<f>}] (6) 
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The Euclidean effective action S{(f)}, corresponding to the Hamiltonian of Eq.(^ has the 
form (for simphcity we ignore charge and magnetic frustration for the time being) 



S[<P] = rfr ^ ^(0.)2 + ^ ^(0, - 0,)2 -EjY: cos(</>, - </>,) (7) 

° . 86 ^.^.^ ^.^.^ J 

with (3 being the inverse temperature. The first two terms are easily recognized as charging 
energies expressed in terms of voltages (see Eq.(^)). In the presence of dissipative tunneling, 
the effective action has a Caldeira-Leggett form and acquires an additional term 

Sd[4>] = \f drdr' ^ a{T - t')F{<P,,{t) - 0,, (r')) (8) 

^ <ij> 

where (pij = (pi — (pj. Both the dissipative kernel ^(t) related to the I — V characteristic 
of tunnel junctions and the function F{{(f)}) depend on the nature of the dissipation. 
Various mechanisms: tunneling of quasi-particles and/or the flow of Ohmic currents through 
the substrate or between the junctions themselves PD| , P5| . 

For an Ohmic junction, as is the case when the bath is formed by quasi-particle excita- 
tions in normal metals (or gapless superconductors), the kernel is 

TT 1 1 

air) 



2emN f3^ sm^nr/p) 



Here, the normal-state resistance Rn controls the coupling to the environment |^9[- For ideal 
superconducting islands on the other hand, the BCS gap inhibits leakage currents at small 
voltages and as a consequence the dissipative kernel is short range in imaginary time. There- 



fore, quasi-particle tunneling results in a renormalization of the junction capacitance PS 



The dissipation mechanism also affects the form of the function F(0) in Eq. (^. If 
normal-electron tunneling occurs via discrete charge transfer, as it is for a quasi-particle 
current, -F(0) is a periodic function of the phase 

FQp(0,,(r)-0,,(r')) = l-cos(^^^^^^Mlj . (IQ) 

If, on the contrary, dissipation is due to normal shunts or more generally to the interaction 
with the electromagnetic environment, F is quadratic in 

F„(^,M-^,M) = i(Ml)^Mi:))' (u) 

indicating that the charge at a junction can assume continuous values. 

One can also consider dissipation due to currents flowing to the substrate. This local- 
damping model plays an important role in classical, proximity coupled Josephson arrays. 
Voltage fluctuations compared to the ground instead of voltage differences between junctions 
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are the crucial variables: the dissipative part of the action now depends on the phase 0j of 
each island and not on the phase difference (pij 

fM^.(r)-A(r')) = i(*M^)' . (12) 

The path integration in Eq.(P) depends on the nature of dissipation mechanism and it is 
related to the charge on the islands being a continuous or discrete variable [jlO[- In the first 
case, the phase is considered an extended variable and in the path integral 0j(O) = 4>i{P)- If 
the charge is a discrete variable, a summation over winding number is implied in Eq.(^) 

D<P — I[ d<P,o E / D<P^i^)- (13) 

i ■'^ {m,=0,±l,...}-^<^'0 

These non-trivial boundary conditions express the fact that the charges of the grains are 
integer multiples of 2e. 

The dissipative coupling strength is usually expressed in the form a = Re/Rq- The exact 
value of the effective resistance R^ is not a priori clear. Consider an array of unshunted 
tunnel junctions. If thermally excited quasi-particles were the only source of damping, a 
measure of R^ would be the subgap resistance which is many orders of magnitude larger 
than the normal-state resistance R^. However, measurements hint at a much smaller Rf. 
that is closer to Rn- The exact mechanism is not clear but one should always keep in mind 
that Josephson junctions are highly nonlinear elements. Some coupling to higher energy 
scales may therefore occur and Re is smaller than the subgap resistance but not lower than 
Rn- This coupling may for instance occur when a vortex crosses a single junction thereby 
producing voltage spikes. Throughout this review, we keep the notation simple and use a 
to denote the dissipation strength regardless of the underlying dissipation mechanism. Its 
origin will be specified from case to case. 



3. Related Models 

The S-I transition has been investigated by studying model Hamiltonians, the so called 
Bose-Hubbard and XXZ models, closely related to the QPM of Eq.(^. We follow the 
notation used currently in the literature and point out the connection with the couplings 
used to defined the QPM. The Bose-Hubbard (BH) model 0] is defined as 



^ = ^ E ^^U,,n, - E ^« - I T.i^~'^'' bib, + h.c.) (14) 

Here, b\ b are the creation and annihilation operators for bosons and rij = bjbi is the number 
operator. Uij describes the Coulomb interaction between bosons {Uij — > E^j'), is the 
chemical potential, and t the hopping matrix element. The connection between the Bose- 
Hubbard model and the QPM is easily seen by writing the field 6j in terms of its amplitude 
and phase and by subsequently approximating the amplitude by its average. This procedure 
leads to the identification bi ~ e**^' . The hopping term is then associated with the Josephson 
tunneling {< n > t — > Ej) while the chemical potential plays the same role as the external 
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charge in the QPM (/x — > Qx)- The mapping becomes more accurate as the average number 
of bosons per sites increases. 

In the case of strong on-site Coulomb interaction Uu oo and very low temperatures 
only few charge states are important. If the gate voltage is tuned close to a degeneracy 
point, the relevant physics is captured by considering only two charge states for each island, 



and the QPM is equivalent to an anisotropic XXZ spin-1/2 Heisenberg model pT] , ^ 



Hs = -hJ2S^ + Y: St U,, S;-EjY: (e^^- St Sr + e-^- 5+ S-) . (15) 

The operators S'f, S^, S~ are the spin-1/2 operators, S'f being related to the charge on 
each island (g^ = + |), and the raising and lowering Sf operators corresponding to the 
"creation" and "annihilation" operators e^""^^ of the QPM. The "external" field h is related 
to the external charge hj h = {q^ — l/2)J2j Uij. Different magnetic ordered phases of the 
XXZ Hamiltonian correspond to the different phases in the QPM. Long-range order in {S~^) 
indicates superconductivity in the QPM while long-range order in (S^) describes order in 
the charge configuration (Mott insulator). 

The three models are equivalent in the sense that they belong to the same universality 
class (they lead to the same Ginzburg-Landau effective free energy). The non- universal 
features like the location of the phase transitions, the shape of the phase diagram (and 
sometimes the very existence of some intermediate phase) depend quantitatively on the 
specific choice of the model. A more rigorous discussion of the different dynamical algebras 
realized by the three models can be found in Refs. p3H47[ . 



A generalization of these models to the case in which amplitude fluctuations are coupled 
to phase fluctuations was discussed in Ref. Very recetly Yurkevich and Lerner |^ 



developed a nonlinear sigma model description of granular superconductors that reduces to 
the BH model in the limit in which amplitude fluctuations are ignored. 



B. The zero-field phase diagram 



Superconductor-Insulator transition has been investigated in great details in J JAs [p0|-[53 

The flrst con- 



as well as in granular systems [54-60| and uniform ultra-thin fllms [31-63 



trolled measurements on the S-I transition in junction arrays have been carried out by 
Geerligs et al. Part of the data of Ref pOl together with the new data of Ref. 151 



j]. Part of the data of Ref 
are presented in Fig. ^. It shows the resistive behavior of six different square arrays in zero 
magnetic fleld. The zero-bias resistance per junction Rq{T) has been measured with a very 
small transport bias (current per junction smaller than < in the linear part of the 

current- volt age characteristics. Three arrays become superconducting; two arrays insulating 
and one array that lies very close to the S-I transition shows a doubly reentrant dependence. 
The horizontal dashed line in Fig. ^ is the critical resistance value of SRq/ti (see Eq. (^)). 

For the three arrays that become superconducting, the data are fltted to the predicted 
BKT square-root cusp dependence on temperature. 



Ro{T)/Rn = cexp{-b[Ej/{T - Tj)] 



1/2N 



with b and c constants of order one. In order to compensate for the temperature dependence 
of Ej, it is convenient to deflne a normalized temperature is deflned as T/Ej{T). From the 
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fits the normalized BKT transition temperature Tj is determined. Near the S-I transition 
Tj is substantially smaller than the classical value of COO-Ej. Note that at low resistance 
levels {Rq{T) < IO'^-Rat), deviations from the square-root cusp dependence are found and 
that the resistance decreases exponentially. This is indicative of thermal activation of single 
vortices across the whole array width |^ . 



The two arrays with a ratio of Ej / Ec < 0.55 show a continuous increase of the resistance 
as the temperature is lowered, i.e., the arrays become insulating at zero temperature despite 
the fact that each island is still superconducting with a well developed BCS gap! It has been 
proposed that, due to the long range interaction between the charges, the conductance will 
follow a square-root cusp dependence on temperature in a similar way as the resistance for 
the superconducting samples. This square-root cusp dependence characteristic for a charge- 
BKT transition is generally not observed. Instead the conductance decreases exponentially 
as temperature is lowered. This issue will be discussed in some details in the section devoted 
to the BKT transition (section II. B. 4). 

The resistance of sample with Ec/ Ej = 1.7 has a very remarkable dependence on tem- 
perature. Starting at high temperatures, Rq{T) first decreases when the temperature is 
lowered. Below T=150mK, however, R{T) increases by more than three orders of magni- 
tude and at the same time a charging gap develops in the I-V curve. Finally at 40 mK, 
Ro{T) starts to decrease again. The second reentrant transition at 40 mK seems to be a 
more general feature of arrays near the S-I transition which is also present in a magnetic 
field. 

Reentrant behavior in the resistance has also been observed in granular superconductors. 
Already in the early theoretical works on the QPM various explanations have been proposed. 
We will summarize some of the ideas. Efetov |^5[ suggested that the thermally excited quasi- 
particles could screen the Coulomb energy thereby lowering the threshold for the onset of 
phase coherence. Stimulated by Efetov's work a number of theoretical papers showed that 
a reentrant phase boundary can be obtained in the QPM however it turned out that the 
very existence of the re-entrance was sensitive to the approximation scheme used. Moreover 
it was shown that even if there is no re-entrance in the phase diagram, the QPM leads 
to a fiuctuation dominated region p6|-p8| which may account for the observed re-entrance. 
Many physical ingredients not contained in the QPM (random offset, dissipation,...) may 
play a role as well. Recently Feigelman et al. proposed that the parity effects may be 
responsible. At intermediate temperatures the screening of quasi-particles would decrease 
the effective Coulomb interaction (and therefore the resistance). At lower temperatures 
screening disappears due to the excess free energy associated with odd grains, leading to an 
upturn of the resistance curve. 

The data of Fig.^ can be used to construct a phase diagram for phase transition of Joseph- 
son arrays in zero-magnetic field |5l[] as shown in Fig. ^. In this figure the superconducting- 
to-normal phase boundary is the vortex-BKT phase transition. Temperature on the vertical 
axis in this figure is given in units of Ej and scaled to the classical (in absence of charging 
effects) BKT transition Tf\ The experimental value Tf^ = 0.95Ej is close to the theo- 
retical value of 0.90 determined from Monte Carlo simulations. On the insulating side of 
the figure no strict phase transition was observed. The dashed line therefore is somewhat 
arbitrary. It represents the crossover to the low-temperature region with Rq > IO^Rn- Fig.|^ 
indicates that at zero temperature the S-I transition takes place at Ec/Ej ^1.7. 
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The existence of a zero temperature (quantum) phase transition can be understood by 
simple arguments as already discussed in the previous section. Mean-field |70], variational 



approaches [[ri|-|73|, 1/z expansion |]73[, Monte Carlo simulations fT^JT^ and cluster expan- 
sions were also applied to the QPM model of Eq. (^. We refer the reader to the various 
articles for a discussion and comparison. Here, we present a few approaches that give a 
self-contained description of the phase diagram. We start with the mean-field calculation. 
Inaccurate in determining the critical behavior, this approach is yet capable to capture 
most of the features measured in the experiments. We first consider the case in which both 
charge and magnetic frustration are absent. Note that for the time being, dissipation is not 
included. It will be considered in detail in Section IlIE. 



1. Mean-Field approach 



The mean-field decoupling consists in approximating the Hamiltonian of Eq.(0) by |65|j70 



Hmf = ^YlQi^i/Qi - zEj(cos0)mf XI • 

i,j j 

The average (cos(0))mf plays the role of the order parameter and it should be calculated 
self-consistently 

(cos(0))mf = Tr {cos(0i) exp(-/3ifMF)} /Tr {exp(-/3ifMF)} • 

Close to the transition point the thermal average on the r.h.s can be evaluated by expanding 
it in powers of (cos(0))mf (This can be done because the transition is continuous). The 
transition line is determined by the equation 

1-zEj / dr {cos (f)i{T) cos (f)i{0))ch =0 (16) 
Jo 

where the average < ... >ch is calculated using the eigenstates of the charging part of the 
Hamiltonian only. Note that the (imaginary) time evolution of the phase is due to charging 
as well. 

In the classical limit the phase correlator is unity (there are no quantum fiuctuations) and 
the mean-field transition temperature is zEj/2. Charging effects inhibit phase fluctuations 
and the critical temperature decreases. Explicit formulas for the phase-phase correlator are 
given in Appendix y. At T = in the self-charging limit {E^^^^ = E^dij) the correlator reads 

(cos0i(r) cos0i(O))ch = ^exp{-4rEo(l -r//3)} . 



By substituting this expression in Eq.(|T^) one gets a value for the S-I transition which 
coincides with the simple estimate based on energy considerations. 

The detailed structure of the phase diagram depends, even in the absence of (magnetic or 
electric) external frustration, on the range of the electrostatic energy. The phase diagram in 
the short-range case is sketched in Fig.^. Very recently a detailed analysis of the dependence 
of the phase boundary on the form of the capacitance matrix has been performed in Ref. 
||78| using perturbation theory and numerical simulations. 
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For two-dimensional arrays the transition to the superconducting state is of the BKT 
type with no spontaneous symmetry breaking. Quantum fluctuations renormalize the value 
of the transition temperature but do not change the universality class of the transition. The 
corrections to the BKT transition due to quantum fluctuations have been evaluated in a 
semi-classical approximation in Ref. ||79[. For Ej » Ec, the JJA behaves as a classical XY 
model but with a renormalized Ej. This approach breaks down when quantum fluctuation 
drive the transition to zero. At zero temperature there is a dimensional crossover and the 
S-I transition belongs to the {d + 1)-XY universality class. 



2. Coarse- Graining approach 

Although it is very useful in determining the structure of the phase diagram, the mean- 
field approach has various shortcomings. For a more accurate description of the quantum 
critical regime one has to resort to different approaches. Universality implies that the critical 
behavior of the system depends only on its dimensionality and on the symmetry which is 
broken in the ordered phase. Many properties of JJAs can be extended to other systems 
that show a S-I transition. By using the coarse-graining approximation it is possible to go 
from the QPM model to a Ginzburg-Landau model with an effective free energy which is a 
function of the order parameter |S^JST|] only. Since the transition is governed by quantum 
fluctuations, the order parameter depends both on space and (imaginary)-time [|T3| , p2| . 

The coarse-graining proceeds in two steps: 

• An auxiliary field ijj{x,T) (which has the role of the order parameter) is introduced 
through a Hubbard- Stratonovich transformation. The partition function is then ex- 
pressed path integral over ^p. 

• The assumption that the order parameter is small close to the transition allows for 
a cumulant expansion to obtain the appropriate Ginzburg-Landau free energy. The 
coefficients depend on the details of the microscopic model. 

The partition function of the QPM is given by 

Z = Tr{e-^^QPM} = Z,,(T,e-/o''^-^^M),, (17) 

where the subscripts ch and J refer to the charging and Josephson part of the Hamiltonian 
in Eq.(^. By applying the Hubbard- Stratonovich transformation to the Josephson term one 
gets (in absence of magnetic and charge frustration) 



exp J — dry e'^^e-'^' + h.c. 

2 JO 7~t 



j Orn^Pexp I - j^''rfr^[i?j]-.V;(r)^,(r) + /^""^rj: [^r(r)e*(^) + h.c] | (18) 



Here we introduced a matrix [Ej]ij which is equal to Ej if i and j are nearest neighbors and 
zero otherwise. The partition function can be written as 
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Z = Z,hj Dtlj*Dt/jexp{-F 



(19) 



Close to the phase transition one can perform a gradient expansion: 

F = J:J drdr'Ur) [\Ej]-15{t - r') - g{r - r')5.,] r,{r') + ^T. J • (20) 

The dynamics of the field ip is governed by the phase-phase correlator 



9[T 



(exp[0j(r) - 0i(r'))c/, 



that was already encountered in Eq. ( [I6| ) . The coefficient k is related to the four-point phase 
correlator. 

The partition function in Eq.(]TP|) can be calculated using a mean-field approximation for 
the phase by evaluating it in the saddle point approximation. The results coincide with that 
of the previous subsection. In the coarse-graining approach, however, a systematic treatment 
of the fluctuations is possible and it allows to study transport as well. In the case of zero 
offset charges and zero external magnetic field, by expanding the phase correlator around 
the zero-frequency and zero-momentum value, the quadratic part of Eq. (pOD can then be 
rewritten as 



where, using the expression given in Appendix |C|, the coefficients can be expressed as: 



7 

c 



g-\ujn = 0) Ej 
g'^jujn = 0) 

4En 



(21) 



At T = 0, this system belongs to the same universality class as the {d+ 1) XY model. One 
can readily obtain all the critical exponents from what is know from the XY model p3[. The 
dynamical critical exponent is z = 1 due to the symmetry between space and time. At finite 
temperatures the transition belongs to the Berezinskii-Kosterlitz-Thouless universality class 
and there is no spontaneous breaking of the symmetry. The dynamical critical exponent 
and the dimensional crossover is modified in the case of 1/r-Coulomb interaction between 
charges [El . 



3. Duality transformations 

Duality transformations have proven to be a powerful tool in field theory and statistical 
mechanics [^|. The idea behind it is that the weak coupling region of a particular system can 



be mapped onto the strong coupling range (and vice versa). The symmetries of the system 
under this transformation lead to important insight into the structure of the model, especially 
in the region of intermediate couplings which is usually elusive to standard treatments. Dual 
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transformations constitute a powerful approach since it is possible to recast the partition 
function solely in terms of the topological excitations of the system [^ -89t . 



In this section we derive some properties of quantum JJAs derived from duality. There is 



a dual transformation |]90|-p2| relating the classical vortex limit, Ej ^ Ec, to the opposite 
charge limit, Ej ^ Ec- The situation is most transparent in the case Cq <^ C, which might 
be more relevant for arrays. The interaction between charges on islands is then logarithmic. 



analogous to vortex interactions in classical, superconducting arrays. The charges form a 
2D Coulomb gas and are expected to undergo a BKT transition at T^^^-* ~ Ec ||9^ (see also 
next subsection). 

Using the results discussed in Appendix the partition function of a JJA can be ex- 
pressed as a sum over charge q and vortex v configurations 

Z = Y, e-^^'^'"^ . (22) 

The effective action S{q, v} reads 

S{q,v} = f drY, [2e\,{T)C;^q,{T)+iTEjv,{T)G,,v,{r) 

+ i qi{r)QijVj{T) + j^qi{T)Gijqj{T)^ . (23) 

This action describes two coupled Coulomb gases. We have used a continuous time notation 
for clarity. Since q and v are integer valued fields, the path integral is well defined on a 
discretized time expression. The charges interact via the inverse capacitance matrix (first 
term). The interaction among the vortices (second term) is described by the kernel Gij, 
which is the Fourier transform of k^"^. At large distances Vij 3> a between the sites i and j 
it depends logarithmically on the distance: 

The third term describes the coupling between the topological excitations in the two limits, 
i.e., it describes the coupling between charges and vortices. The function 

. ( Vi-yA 

Kyij = arctan | | 



represents the vortex-phase configuration at site i when its center is placed at the site j. 
The coupling has a simple physical interpretation: a change of vorticity at site j produces a 
voltage at site i which is felt by the charge at this location. The last term qCq stems from 
the spin-wave contribution to the charge-correlation function. In the limit in which Ej — »• 
or Ec the action for a classical system of Cooper pair charges or vortices is recovered. 

The effective action in Eq.(p3D shows a high degree of symmetry between the vortex and 
charge degrees of freedom. In particular, in the limit Cq <^ C the inverse capacitance matrix 
has the same functional form as the kernel describing vortex interactions: 

e Oj- — '^ij ■ 
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Hence charges and vortices are dual. There is a critical point for which the system is self-dual 
with respect to interchanging them: 

Ec vr2 

The duality is strict for vanishing self-capacitance and in the absence of the spin-wave duality 
breaking term (qGq) in Eg . (p3D . This latter term is irrelevant at the critical point, i.e., it 
merely shifts the transition point. However, it has important implications for the dynamical 
behavior. 

Duality transformations have also been applied to a three-dimensional JJA consisting 
of two 2D arrays placed on top of each other |p^-p6|. The authors of these papers assume 
that there is only capacitive coupling between them (no Josephson coupling). The most 
interesting situation arises when one array is in the quasi-classical (vortex) regime while the 
other is in the quantum, charge regime. Then, vortices in one layer and charges in the other 
one are well-defined. 

The Euclidean effective action, in term of phases, ^2} has the form 



'^{01,02} = dr 



8e' 



''^i^i - 4>jii)^ - Ejf, cos(0i^ - (f)jf,) 



(24) 



where Cq/^ are the island capacitances in array fi relative to ground, are the junction 
capacitances in the array fi, and are the interlayer capacitances between islands on top 
of each other, while Ej^ are the Josephson coupling constants in the layers. Similar as in a 
single array, we move from a description in terms of phases to one in terms of charges and 
vortices, and use the duality of the resulting action to investigate the transition. 

Before we proceed with the calculation it is necessary to stress that in the regime of 
interest the interlayer capacitances not only couple the layers, but also renormalize the 
island capacitances Cqi and C02 to ground. The physical reason for this is that due to the 
strong fluctuations of charges in layer 2 and vorticities in layer 1 these variables are effectively 
continuous, and hence a coupling to the other array plays the same role as a coupling to 
ground. Due to screening, the interaction between charges in each layer has a finite range 
for any non-zero C^, and the BKT transition is replaced by a crossover. However, in the 
limit Cqi ^ Cx ^ Ci the screening length ,^1 ~ (Ci/C^Y^'^ can be large enough to make it 
meaningful to speak about the charge-unbinding transition (the transition is exponentially 
sharp). Below we consider this case. For not so weak coupling, on the other hand, this 
description becomes meaningless, since the crossover is strongly smeared, and the insulating 
phase is absent. 

It is possible to introduce vortex degrees of freedom in the same way as for one array. 
We obtain the effective action for charges qair) in layer 1 and vorticities Vi2{T) in layer 2 
(to be referred below as and Vi) 



S{q,v}= [ dr I J2 (li{'^)Gijqj{T) + , \ ^i('^)Gijqj{r) + ■kEj2 ^ Vi{T)G 

Jo I TC ■■ 47ri!/7i 
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»J^C2 ij L 4vr C1C2 J ivrCi .^.^ J 

This equation looks rather similar to the effective charge-vortex action in one Josephson 
junction array (see Eg. p3|) ) . The most important difference is that while in one layer either 
charges or vortices are well-defined degrees of freedom, Eq. (^5]) describes the system of two 
well-defined dynamic variables on each site - the charges in layer 1 and the vortices in layer 
2. The action shows a duality between charges and vortices (the second term in the square 
brackets is small for <^ Ci,C2). Both kinetic terms for charges and vortices violate the 
duality due to the numerical coefficients. However, close enough to the transitions these 
terms only produce a small renormalization of the transition temperature, and are therefore 
irrelevant. Another interesting feature of this action is that the last term, describing the 
interaction between charges and vortices, is also small, while in a single-layer array the 
interaction is always of the same order of magnitude as the other terms. 



4- Berezinskii-Kosterlitz-Thouless transitions 

In classical arrays, it is well established that JJAs undergo a BKT phase transition 
to the superconducting state. In the opposite hmit where charging dominates, quantum 
fluctuations of the phases are strong, and vortices are ill-defined objects. In this regime the 
duality transformations discussed in the previous section show that charges on the islands are 
the relevant variables. Similarly to vortices in classical arrays, they interact logarithmically 
with each other and are expected to undergo a charge-BKT transition leading to insulating 
behavior [pSj . 

A critical point separates the superconducting and insulating regime at T = 0. As 
discussed before, the various models give different estimates for the value of the critical 
point, but in all cases it lies close to Ej/Ec ~ 1. The theoretical phase diagram in the limit 
of logarithmically interacting charges is shown in Fig.^ For any T ^ the array has three 
phases. Next to the superconducting and insulating phase, there is a region with normal 
conduction. Here, there are always some free vortices/ charges present as they are generated 
by thermal or quantum fluctuations. Experimental verification of this diagram has been 
reported in Refs. [0,0 In the remainder of this subsection, we discuss the experimental 



aspects of BKT transitions in arrays. 

The BKT transition in the classical case has been studied in great detail (see e.g. 
Refs. On approaching the critical point of the S-I transition, the vortex-BKT transi- 

tion temperature is lowered by quantum fluctuations. For Cq = the shift of the transition 
temperature is 



r.-^fl-^#) . (26) 



37cE 



A reduction of the transition temperature is generally observed in quantum arrays. The 
Delft data (shown in Fig. however, show that with increasing Ec/Ej ratio the reduction 
of the transition temperature goes faster than predicted above (Eq. 
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On the charging side {Ej <^ Ec) of the phase diagram, the charge-BKT transition occurs 



at a temperature |^ 



rp Ec 
J-rh — 



TT 



E] 
0.31-^ . 



Note that the charge-BKT transition exists for both arrays with superconducting and normal 
islands ||100|| . 

The existence of a charge-BKT transition implies that arrays are insulating below the 
transition temperature. The conductance should vanish with a characteristic square-root 
dependence according to 



G ~ RJ,^ exp {-26 [T/T,^ - 1]"^/'} 



(27) 



where the constant 6 ~ 1. The temperature dependence of the conductance in the charge 
regime has been investigated by several groups [p.01H103|| . Instead of the predicted square- 
root cusp behavior, an exponential (activated) temperature dependence 

has been observed with an activation energy 



A + QMEc 



In arrays the screening length is about 10^ lattice constants and therefore there is no 
a priori reason for observing such dramatic deviations from the BKT theory. Recently, 



Feigelman et al. |^9[ re-examined the problem and found that parity effects together with the 
screening of the Coulomb interaction due to thermally activated quasi-particles is responsible 
for masking the charge-BKT transition. At temperatures above a crossover temperature T* 



where parity effects ||104|| set in, the presence of quasi-particles rules out the possibility to 
observe the charge-BKT transition. If the BKT transition temperature is larger than the 
crossover temperature, the array behaves as a normal one and a charge BKT transition occurs 
associated with the unbinding of quasi-particles at temperatures close Ec/4:7i. The presence 
of free charges screens the interaction between Cooper pairs resulting in the unbinding of 
pairs. The resistance as calculated by Feigelman et al. is expected to be 



In 



R{T) 



R{Ec/47r] 



~ mm 



Fp{T) 



E, 



c 



1/2 



where Fp{T) is the free energy difference between islands with an even and an odd number 
of electrons [|104|| . From analogous considerations one may conclude [39| that normal arrays 
are better suited for studying the charge unbinding transition. 



C. Magnetic frustration 

When in the classical limit {Ec << Ej) a perpendicular magnetic field is applied, vortices 
enter the array above a certain threshold field |^4|. Just as in films, the vortex density 



increases with increasing magnetic field. In junction arrays the periodic lattice potential, 
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however, prevents vortices to move at low temperatures. Only above the depinning current, 
vortices move (flux- flow branch). The resistance in this branch increases approximately 
linearly with / = $/$o up to / ~ 0.2 ($ is the magnetic flux piercing through an elementary 
plaquette). A phenomenological model, analogous to the Bardeen- Stephen model used to 
describe flux-flow in fllms, is in good agreement with experiments providing that coupling 
between spin waves and vortices is taken into account (see section lll.C) The properties of 
arrays at low frustration (low vortex densities) are dominated by single-vortex properties 
and are the subject of the next chapter. 

In larger magnetic flelds commensurability effects come into play and the behavior of 
junction arrays is richer than that of fllms. A magnetic fleld applied perpendicularly to the 



array leads to frustration ||106|] . The presence of the magnetic fleld induces vortices in the 
system and if the frustration is a rational number, f = p/q, the ground state consists of a 
checkerboard conflguration of vortices with aqxq elementary cell. The stability of the vortex 
lattice against a bias current leads to a decrease in the small-bias resistance at flnite fractional 
flllings. In order of their relative strength, one expects dips at / = 1/2, 1/3, 1/4,2/5, ... in 
square arrays and at / = 1/2, 1/4, 1/3, 3/8, ... in triangular arrays as is illustrated in Fig.||. 
Near these fractional values of /, defects from the ordered lattice (excess single vortices or 
domain walls) are believed to determine the array dynamics in a similar way as the fleld 
induced vortices determine array dynamics near / = 0. Therefore, arrays near commensurate 
values with high stability such as / = 1/2 may qualitatively behave in a similar way as near 
zero magnetic fleld. Because all properties are periodic in / with period / = 1 an increase 
beyond / = 1/2 does not lead to new physics. 

A particularly interesting case is the fully frustrated situation (/ = 1/2) in square arrays. 
The two degenerate ground states consist of a vortex lattice with a 2 x 2 elementary cell. The 
current corresponding to this vortex arrangement flows either clockwise or anti-clockwise in 
each plaquette (chiral ground state). Interaction between domain wall excitations with 1/4 
fractionally charged vortices ( at the corners of a domain wall) and excess single integer 
vortices are believed to trigger a combined BKT-Ising transition. A fully frustrated array 
has two critical temperatures related to the Z2 and U{1) symmetries of the problem. Their 
existence has been investigated both by analytical methods and Monte Carlo simulations. 
Even at the classical level, the complete scenario is not fully understood yet. There is 
numerical evidence either supporting the existence of two very close critical temperatures 
with critical behavior typical of Ising and BKT transitions respectively or the existence of a 
single transition with novel critical behavior. Further reference to classical frustrated arrays 
can be found in Ref. [|107]| . 

At the mean-fleld level the full phase-diagram including charging effects and magnetic 
flelds is obtained by solving an eigenvalue problem equivalent to the Hofstadter prob- 



lem ||108|| . The resulting phase boundary as a function of vorticity shows commensurability 
effects. Although the superconducting transition temperature is reduced, the average conflg- 
uration of the phases and the supercurrent flow patterns are unchanged. The ground state is 
still chiral |p,09|| . More detailed calculations based on expansion in Ej/Eq of the QPM | 11U ] 
and on the BH model ||111|| conflrm the butterfly-like behavior of the S-I transition. In Fig.^ 
the theoretical results obtained by Kim et al. [|11(]|| are shown. 

As in the unfrustrated case, measurements indicate a S-I transition at T = 0. For the 
same set of samples as presented in Fig.^ the S-I transition for / = 1/2 has been studied. 



22 



The phase diagram is shown in Fig.|TO|. The transition takes place very close to a normal- 
state resistance of 11 kf2. The critical Ec/Ej ratio is about 1.2, a factor 0.7 lower than 
the zero-field value. This decrease of is consistent with the simple model that involves a 
reduction of effective Josephson coupling energy at / = 1/2: the interaction energy of a 
vortex pair is a factor \/2 smaller than in zero field. With this lower effective coupling the 
critical value of Ej / Ec of the S-I transition is reduced by a factor y^, which is close to the 
observed reduction of 0.7. 



The experimental data agree rather well with the quantum Monte Carlo calculations |112 



in the classical limit. The experimental points of the transition temperatures are, however, 
lower than the calculated ones by entering in the quantum regime. At present, there is 
no explanation for this discrepancy. It would require a more detailed study and better 
understanding of the phase transition at / = 1/2. The calculations do indicate, on the other 
hand, a S-I transition at Ej/Ec ~ 1, in agreement with the experiment. 



In Fig. 11, the critical Ec/Ej ratio as a function of applied magnetic field for square 
arrays is plotted. The three points at / = 0, 1/2, and 1/3 are combined with two data 
points of the field-tuned S-I transition (see Section [11 H[ ). After a rapid decrease the critical 
ratio is almost constant for / > 0.1. The critical Ec/Ej ratio at / = 1/2 is larger than at 
other nearby values of /, indicating once again the stability of the phase configuration at 
/ = 1/2. Fig. 11 also indicates that arrays in the range 1.2 < Eq/Ej < 1.7 do not show 
special behavior at commensurate /-values (e.g. dips in the magnetoresistance) ; arrays are 
superconducting in zero field but insulating at / = 1/2, 1/3, 1/4, ... 

It is possible to derive a Ginzburg-Landau free energy also in the presence of rational 
frustration ||113|| . The calculation proceeds along the lines outlined for the / = case in 
Section II B 2 . The important difference is that one should expand the free energy about the 
most fluctuating modes. In the unfrustrated case this means an expansion about k = (0, 0). 
For the fully frustrated case (/ = 1/2) the expansion is carried around the two points: 
k = (0, 0) and k = (0, vr) thereby reflecting the superlattice structure of the ground state. 
The resulting free energy depends on a multicomponent (complex) order parameter (e.g. 
two coupled complex fields in the fully frustrated case). 

The magnetic-field dependence of the critical exponent zu {u governs the divergence of 
the correlation length and z is the dynamical critical exponent) was considered by Niemeyer 



et al. ||111|| . In zero field the mapping onto a three-dimensional XY model implies that 
zv = 0.67. Their analysis hints to a dynamical exponent that increases with the magnetic 
field. It is, however, difficult draw conclusions on the values of zu. It could smoothly increase 
with magnetic field, or immediately jump to one on once the magnetic field has switched 
on. Combining the fact that zu < 1 for / = 1/2 and that a higher-order expansion predicts 
zu = 1 there, the authors conclude that the answer is zu = 1 for all nonzero magnetic fields. 
Monte Carlo simulations by Cha and Girvin [|1 14|| obtain for the / = 1/2 and / = 1/3 cases 
the values z = 1 and 1/z/ = 1.5, consistent with the analysis outlined above. 

Finally we mention that, in addition to two-dimensional arrays, frustration effects can 



be studied in quantum ladders as proposed in Refs. ||115 
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D. Charge frustration and the supersoHd 



A uniform charge can be introduced in a quantum JJA by applying a gate voltage 
with respect to the ground plane. This effect is known as charge frustration. Although 
from a theoretical point of view charge and magnetic frustration are dual to each other, 
experimentally it is only possible to tune the magnetic frustration in a controlled way. In 
all arrays random offset charges, presumably caused by defects in the junctions or in the 



substrate ||116|| , are present. Electron or quasi-particle tunneling will partly compensate 
these offset charges so that their value lies between — e/2 and +e/2. These charges can, 
in principle, be nulled out by the use of a gate for each island; this procedure, however, 
works only for small networks. In large arrays they cannot be compensated because too 
many gate electrodes would be necessary, requiring too complicated fabrication procedures. 
A uniform charge frustration has therefore not been realized yet in 2D Josephson arrays. 
Lafarge et al. [ 117| ] have investigated charge frustration by placing a gate underneath a 



Josephson array. They managed to obtain a 40% variation of the resistance between the 
unfrustrated and the (nominally) fully frustrated array. But most importantly, in studying 
the cur rent- volt age characteristics, it was impossible to quench the Coulomb blockade as 
it can be done in circuits with few junctions. In future arrays, charge frustration may be 
applied more uniformly if the influence of offset charges can be drastically lowered. 



1. Phase diagram 

The energy difference of two charge states in each island with q and q + 1 extra Cooper 
pairs may be reduced by changing (which means to change the external charge q^). Con- 
sequently the effects of a finite charging energy are weakened and the superconducting region 
in the phase diagram is enlarged. It turns out that for certain values of the gate voltages the 
energy difference vanishes implying that the Mott gap, and therefore the insulating behavior, 
is completely frustrated. At the degeneracy point even a small Josephson coupling makes 
the system superconducting since there is no pay in energy for moving a charge through the 
whole array. In general it is intuitive to expect that the extension of the insulating lobe will 
be maximum at integer values of the external charge since in this case the excitation energy 
is highest. 

A quantitative analysis of this phenomenon has only been obtained in models with a 
short-range electrostatic potential. Uniform charge frustration gives rise to two new effects: 

• lobe structures appear in the phase diagrams 

• new states in the phase diagram may occur (Wigner-like crystals and the Supersolid). 

The remainder of this subsection is devoted to the lobe structures and the Wigner-like charge 
ordering; the next subsection treats the supersolid. 

The lobe-structure already follows from a mean-field analysis with on-site interaction 
only. The corresponding phase diagram can be obtained by evaluation of the correlator 



given in Eq.([16D in the presence of an uniform charge. In Fig.|12|, the mean- field phase 
boundary in the presence of a uniform background charge is shown. The detailed structure 
of the lobes is very sensitive to the model used (QPM, BH, XXZ) and on the approximation 
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made [^ , p.l8| - [123|| . The lobes in the Monte Carlo calculation are sharper than predicted 



by mean-field, but smoother than expected from the strong coupling expansion. 

In the case of a finite-range interaction a number of new classical ground states exists 
characterized by a crystal like ordering of charges. The phase diagram contains extra lobes 
where the charge density is pinned to a given fractional filling. Their range of existence, in 
the limit of vanishing Josephson coupling, is found by minimizing the charging part of the 
Hamiltonian for a given charge filling. We illustrate this by considering the simple case of 
short range charging energy (only on-site and next-neighbors). In the case of > > 1/2 
(and for square lattices) only three different charge configuration should be considered. 

• All the islands are neutral 

• A checkerboard state can be formed in which a sublattice is neutral and the other is 
charged with one extra Cooper pair 

• All the islands can be uniformly charged with charge 2e. 

The corresponding energies of the different ground states are respectively 

E,H,ml^N = Eoql + ^E,ql (28) 

Echfiif^N = ^Eoql + ^Eo{l - q.,f - ^-E^q.^l - q.) (29) 
EcHmI^N = Eo(l - q^f + |Ei(1 - q^f (30) 

where is the number of islands in the array. The ground state energy is given by Ech,Qo 
for > > qx,i, Ech,Qi for q^^i > qx > <lx,2 and Ech,oo for larger q^. The critical values at 
which the ground state changes are 

Qx,i 

(lx,2 



1 1 



2 1 + {z/2)Ei/Eo 

1 1 + zEi/Eq 

2 1 + {z/2)EjE^ 



The S-I boundary (for finite Ej) can be determined, for example, using a mean-field 
approach. The result is presented in Fig.[T^. The checkerboard state for q^^i > qx > (lx,2 
can be thought of as a Wigner crystal of Cooper pairs. The role of a longer-range charging 
energy (next nearest neighbors,...) is to stabilize the crystalline phases with lower fillings 
(1/4, 1/8,...). Not all the lobes extend down to the Ej = axis. First-order phase transitions 
between different checkerboard states are then possible. Note, that the presence of charge 
ordering, characterized by a periodicity 27i/kx,27i/ky, is detected by studying the structure 
factor 5* at a give wave-vector 



s{kx,ky) = —{Y.^'''''mo). (31) 



1 

'J 

A uniform charge frustration changes the symmetry properties of the system. At fc^; = 
the energy cost to create (or remove) a Cooper pair in a given island is the same. The system 
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possesses particle-hole symmetry. For generic values of the external charge this symmetry 
is broken (i.e., i?c/i,oo 7^ Ech,ii)- In the phase diagram shown in Fig. 13 the tips of the lobes 
correspond to a particle-hole symmetric case while away from the tips the symmetry is 
broken. This change of symmetry is reflected as a new term in the quadratic part of the 
Ginzburg-Landau free energy. This new contribution is 

XJ2J dTij*{T)drUr) (32) 

i 

where 

^ = ' 2Eo • 

The coefficient A vanishes in the particle- hole symmetric case (see Appendix The 
particle-hole symmetry has important consequences for the critical behavior of the sys- 
tem The dynamical critical exponent z changes from z = 1 a.t the tip of the lobes to 



z = 2 in the generic case. 



2. Supersolid 

A solid phase is characterized by charges being pinned on the islands whereas a superfiuid 
phase is characterized by phase coherence over the whole system (i.e., charges are delocal- 
ized). At the end of 60's ||124H126| , ^ it was suggested that, in addition to the solid and 



superfiuid phases, a new state should appear, characterized by the coexistence of off-diagonal 
(superfiuid) and diagonal (charge-crystalline) long-range order. This phase is known as the 
supersolid. If vacancies in a quantum crystal such as solid ^He Bose-Einstein condense, 
they do not necessarily destroy the crystal structure and they may form a superfiuid solid 
(or supersolid). Experiments have been performed on ^He, but no positive identification 
of this coexistence phase has yet been made. There are, however, hints that such a phase 
exists [ p7|JT2l . 



An exciting possibility that attracted a lot of attention was the idea to observe supersolids 
in Josephson arrays [^ , |129| - [133|| . The supersolid phase is located in an intermediate region 
around the half-filling lobe. A simple way to understand its existence is to focus on a region 
close to the phase boundary at ~ 1/2. At densities corresponding to half filling the 
particles form an incompressible solid. Away from half-filling vacancies in the charge-solid 
appear. As they have a bosonic character, they can Bose condense, and therefore they are 
able to move freely through the system. For a small enough density of vacancies one expects 
that the crystal order is not destroyed. 

In the limit of very large on-site charging (hard-core limit in which the island charge 
can only be zero or one) the existence of the supersolid is related to the finite next-nearest 
neighbor interaction as it does not exist for nearest-neighbor interaction only. Furthermore, 



there is no supersolid phase at exactly half-filling. In Fig.|T^ the mean-field phase diagram 
in the hard-core limit is shown. The supersolid region appears in a tiny region away from 
half filling between the superconducting and Mott insulating phases. 

If higher values of charge are allowed, the supersolid phase already exists for nearest- 
neighbor interaction and also at half-filling on the tip of the checkerboard lobe. This is 
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related to excitations which are forbidden in the hard-core hmit. A large nearest-neighbor 
interaction or small on-site interaction favors the supersolid, whereas in the hard-core limit 
the supersolid is suppressed. Thus, it seems that the system itself generates the defects 
(particle-hole excitations) that Bose condense, thereby turning the solid into the supersolid. 
The phase diagram for soft-core bosons is shown in Fig.p!5|. 

Since the supersolid phase is very sensitive to fluctuations, it was important to obtain 
independent checks of its existence. Monte Carlo simulations on the QPM | ]130| | and the BH 
model ||132|| have confirmed the qualitative picture discussed above. In Fig.ll5| the symbols 



represent the phase diagram as obtained from Monte Carlo simulations by van Otterlo and 
Wagenblast ||129| , |130|| . Note that the supersolid region is considerably reduced as compared 
to the mean-field estimates. 

By changing the electrostatic interaction new phases named collinear supersolid were 
found by Scalettar et al. [p.32|| and by Frey and Balents [ p.33|| . A detailed analysis of various 
supersolids including striped phases has recently been performed by Pich and Frey ||134||. Su- 



persolid phases in frustrated systems have been studied as well either by combining the effect 
of charge and magnetic frustration ||135|| or by considering arrays on Kagome lattices ||136| . 

Several other kinds of coexistence phases were studied. The possibility of a spontaneous 
vortex anti-vortex lattice in superfiuid films was explored in Refs. | |137| , |138|| and a coexistence 
phase of superfluidity and hexatic orientational order was proposed in Ref. ||14CI|| . Orienta- 
tional order in incompressible quantum Hall fluids is discussed in Ref. ||141|| . Finally, we 
mention the relation between 2D bosons and 3D flux-lines in type II superconductors (high- 
Tc materials) in a magnetic field ||143| , |144|| . Also in these systems different kinds of long-range 
order may coexist and the equivalent of the supersolid is discussed in Refs. ||145| , |146|| . Re- 
lated is the question whether or not vortices may form a disentangled liquid, which would 
imply a normal ground state for bosons with long-range Coulomb interaction. 



E. Dissipation induced S-I transition 

The behavior of a single Josephson junction with Ohmic dissipation has been discussed 
in a pioneering work by A. Schmid ||147|| who found that there is a zero-temperature phase 
transition governed by the dissipation strength a. Above the critical value a = 1, dissipation 
suppresses quantum fluctuations thereby restoring the classical behavior with a finite critical 
current. For weak damping, on the other hand, quantum fluctuations destroy global phase 
coherence. The supercurrent is suppressed to zero and the junction is in the insulating state. 

Experimentally, this transition has only very recently been detected by Penttila et 
al. | [14(j| ]. The reason is that for a single junction the high-frequency coupling to the en- 
vironment determines the effective damping. Consequently, the effective impedance is of 
the order of 100 [], which is about two order of magnitude smaller than the quantum 
resistance. Penttila et al. increased this impedance, i.e., the decoupled their single junction 
from its environment by placing high-Ohmic, Chromium resistors in the leads close by. 

Dissipation plays an important role in quantum phase transitions of JJAs as well. Orig- 
inally the interest was stimulated by the idea that dissipation could be responsible for 
the observed critical resistance at the S-I transition in arrays and granular films. Later 
Fisher ||150|| and Wen and Zee ||151|| pointed out that the observed critical resistance is a 



zero temperature property associated with the existence of a quantum phase transition and 
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it is not related to the presence of an "extrinsic" source of dissipation. The next section 
discusses the critical behavior on transport properties in more details. Here, we discuss the 
influence of dissipation on the phase diagram. 

The coupling to a dissipative bath has the effect to suppress the quantum fluctuations 
of the phase, i.e. to quench the insulating region. The properties of the environment and 



the type of dissipation are important ingredients in Eq.([lOD and in Eq.(|TT]). As stated 
before, various sources of dissipation should be considered for JJAs. Quasi-particles may 
still play a role at mK as they may be generated by the environment or the motion of 
vortices themselves. From a theoretical point of view and in view of the recent experimental 
advances, it is also possible to realize arrays in which Ohmic shunts are important. These 
shunts can be normal wires placed parallel to the junctions in a a similar way as was done for 
a single junction. Ohmic shunts may also arise because of coupling to the substrate (local 
damping). Very recently the group of Kobayashi succeeded in fabricating a JJA in which 
each junction is shunted by a Cr resistor ||152|| . A different and controlled environment 



was investigated in the experiments of Rimberg et al. [|153|] by placing a 2D electron gas 
underneath the Josephson array. We briefly discuss all these sources of dissipation. 



1. Quasi-particle dissipation 



When the mechanism responsible for dissipation is quasi-particle Tunneling, the effective 
action is that given in Eq.(^). Theoretically this model was studied by means of mean- field 
calculation |p.54| , |155|| , variational approaches ||156| - |159| , and Monte Carlo simulations ||160 



The dominant effect coming from quasi-particle tunneling enters in the renormalization of the 
effective junction capacitance given in Eq.(^. In the mean-field calculation, this amounts to 
a modification of the capacitance matrix in the evaluation of the phase correlator. The zero- 
temperature phase boundary (for short range Coulomb interaction) obtained by Chakravarty 
|154|| , is given by the expression 



et al. 



1 = — a — 
4 In 



1 + ^^a^ 



4 E, 



It is important to stress that in the case of quasi-particle tunneling the array is either 
insulating or superconducting. The interplay between the long-range Coulomb interaction 
and quasi-particle dissipation has been discussed in Ref. [ JT . 

In Section |IIB| , we have interpreted the S-I transition as being driven solely by Coulomb 
interactions. However, given the uncertainty in the damping resistance (e.g. sub-gap resis- 
tance of normal-state resistance) of the junctions the possibility that the transition is driven 
by quasi-particle dissipation cannot be ruled out. The data do not exclude the possibility 
that the S-I transition is influenced by the normal-state resistance. In fact, the Chalmers 
and the group of S. Kobayashi and collaborators have interpreted their 
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group 

data in terms of a Schmid-like diagram. In Fig.|16| we show the results from the Chalmers 
group. The normal state resistance Rn is used as the resistance determining the dissipation 
parameter. If this resistance is used, a reasonable agreement with the theoretical models is 
obtained. 
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2. Ohniic dissipation 



The influence of Olimic sliunts on tlie pliase diagram lias been intensively investigated 
as well. A new phase with local phase coherence is possible, i.e., phase coherence only exists 
as a function of time. Various theoretical methods have been applied in this case such as 
coarse graining [p.62| - |164| , variational ||163| , |165H168|] and renormalization group [|163| , |166| , |169 



approaches. The general trend is, as expected, that the critical value of Ej/Eq for the onset 
of phase coherence is lowered. The dependence on the dissipation strength is stronger as 
compared with the case of quasi-particle damping. As for a single junction, a true dissipative 
transition occurs. A rigorous analysis can be performed in various limits in the Ej — a phase 
diagram | |163| , |166|| . For simplicity we discuss only the T = case and follow the discussion 



presented in Ref. [|I7i 

• In the large a limit time-like fluctuations of the phase are strongly suppressed and 
they only contribute to the renormalization of the effective Josephson coupling. The 
system behaves like a classical JJA with an effective Josephson coupling 



Eyf = Ej(l--\n + - 



At zero temperature the array is in the superconducting phase independent on the 
ratio Ej/Eq. 

• In the case Ej /Eq ^ 1 there is a phase transition at 

1 

a = - 
z 

which separates two phases that both exhibit long-range coherence. Evidence of such 
a phase transition could be detected by measuring the voltage noise power spectrum. 

• In the limit of small damping, the critical ratio of Ej/Eq is renormalized to smaller 
values indicating that dissipation enhances the superfiuid phase. 

• If the ratio Ej/Eq is very small, a dissipative transition to a phase with local order 
can take place at a critical value of dissipation given by 

2 

a = - 
z 

Very recently Takahide et al. [|152|] fabricated a JJA in which each junction was shunted 
by a Cr resistor. By varying the resistance of the shunts and the ratio Ej/Eq they were 
able to map out the phase diagram. The results are in good agreement with the theories of 
dissipation induced quantum phase transition discussed above. 



3. Local damping 



Dissipation may also arise from coupling with the substrate by means of what is named 
as the 'local damping' model. Local damping changes the universality class of the S-I tran- 
sition [p.74|] , and influences the low- frequency dispersion of the vortex response in classical 
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arrays ||1 72| , |1 73|| . As shown in Eg. ([I2D ) , the local damping model correlates the island phase 
in time. In proximity-coupled arrays, which consist of superconducting islands on top of 
a metallic film, the model with local-damping is appropriate to describe the flow of nor- 
mal electrons into the substrate. Aluminum tunnel junction arrays are always placed on 
insulating substrates so that it is not appropriate. 

Dissipation due to local damping is associated with the phase rather with the phase 
difference — as in the resistively shunted junctions (RSJ) model. The number of Cooper 
pairs in each island is allowed to decay, whereas the RSJ model describes only charge transfer 
between neighboring islands. By going through the same steps outlined in the section on 
the coarse-graining method, it is possible to obtain also in this case an effective Ginzburg- 
Landau free energy. The only difference is that now the phase-phase correlator g{T) has to 
be evaluated including the local damping term. For small frequencies the Fourier transform 
reads (for more details see Ref. ||174|| ) 



2 

g{uj^) = g{0) - r]\uj^\' ~ C^l with s = 1. (34) 

a 

The coefficients rj and ( can be determined from the phase correlator (their value is not 
important for our purposes). Using this expression for g{ujn), the free energy in Eq.(pOD 
contains a non-Ohmic dissipative term (oc Icu^H (reducing to Ohmic, or 'velocity propor- 
tional' damping only in the special case s = 1). This means that Ohmic damping in the 
quantum phase model yields non-Ohmic dynamics for the coarse-grained order-parameter. 
The phase boundary in the saddle point approximation is shown in the inset of Fig.p!7|. By 
increasing damping strength, the superconducting region is enlarged. At T = a quantum 
phase transition is ruled out beyond the critical value a = 2. 



4- Tunable dissipative environment 

A controlled study of the dissipative S-I transition has been performed by Rimberg et 
al. ||153|| . They placed a Josephson array on top of a two-dimensional electron gas (2DEG). 
Junction parameters are chosen such that in the absence of the 2DEG the array is insulating. 
The array is capacitively coupled to the electron gas and its screening currents provide 
a source for dissipation. By tuning the back-gate voltage, the electron density and the 
sheet resistance Rg of the 2DEG are varied without changing the array parameters. As the 
resistance of the 2DEG increases the current-voltage characteristics of the array change from 
superconducting to insulating with a Coulomb gap as illustrated in Fig.^. Moreover the 
resistance of the array is very sensitive on Rg as shown in Fig.[l^. Note that in the experiment 
the island capacitance (to the 2DEG) exceeds the estimated junction capacitance of 0.5 fF 
by a factor of 6. 

Wagenblast et al. ||171|| analyzed these measurement and modeled the experimental setup 
by capacitively coupling the array to the 2DEG. Assuming Ohmic dynamics of the 2DEG 
they obtained the following (Caldeira-Leggett like) effective action for the array 

SoM = \llj d\^D,\k,uj,)\^u,uf + Sj (35) 

{Sj is the action related to the Josephson coupling) with the propagator 
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where Co now represents the capacitive couphng to the 2DEG and where I/Qq = RgCo- 
The effective action for the array is Ohmic only in an intermediate frequency range. At the 
lowest and highest frequencies the dynamics is capacitive. The two energy scales are well 
separated in the case Co 2> C and a quantum phase transition is driven by the action at the 
lowest frequencies. As the dissipative action is cut off at the lowest frequencies, a dissipation 
driven transition cannot occur in the strict sense. However, quasi-critical behavior can be 
observed at temperatures and voltages exceeding the low energy scale Qq- In the limit 
fio ^ (Co oo) this behavior converges to a true dissipation-tuned transition. 

An analysis of the conductivity as a function of Ej and a suggests a phase diagram of 
the type represented in Fig.^. The insets show Ro{T) as a function of the temperature 
in different regions. The experiments of Rimberg et al. belong to the right-lower sector 
of the Schmid diagram. The theoretical temperature dependence of the resistivity Rq{T) 
as well as the exponential relation between Rq and Rg are in good agreement with the 
experiments [p.53|| . 



F. Transport Properties 

The unique feature of quantum critical points is that quantum fluctuations, which drive 
the system through the transition, govern its dynamical behavior. The interest in under- 
standing charge transport near a quantum phase transition goes beyond the study of JJAs. 
Important examples are the transition in quantized Hall systems |T^, localization in Si- 
MOSFETs ||175|| and the quantum critical point in cuprates [|176|| . 



In two dimensions right at the S-I transition the zero-temperature conductance ||177 
has been predicted to be finite and universal. This metallic behavior is present even in 
the absence of extrinsic dissipation and it is entirely due to the presence of collective modes 
which become soft at the zero-temperature transition point. Universality at a quantum phase 
transition then implies that the properties of the system are governed by a set of critical 
exponents. X.G. Wen [ [I78|| employed a scaling theory of conserved currents at anisotropic 
critical points identifying universal amplitudes. One of these amplitudes in two dimensions 
reduces to the universal conductance. 



A very simple argument [179| leading to a finite and universal conductance at zero 



temperature can be discussed using the duality between charges and vortices formulated in 
Section p.lB| . Strictly speaking it applies to the case Cq = 0, i.e. for logarithmic interacting 
charges. From the Josephson relation the voltage drop across an array is given by the rate 
of vortices crossing the sample boundary. The current is given by the number of Cooper 
pairs which flow through the system per unit time, i.e., 

y = A(^) I = 2e{q) 

At the self-dual point {v) = (q) and therefore the conductance at the transition (denoted 
with a*) is finite, universal and corresponds to the quantum of resistance for Cooper pairs: 

. 4e2 
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The value of a* changes in the case of short-range charging and/or in the presence of disorder. 
Nevertheless it remains universal (independent on the sample parameters). 

A large amount of theoretical work has been devoted to the determination of the critical 
value a* and the scaling behavior of the conductance. The universal conductance in a 
model with no disorder was considered in Ref. [ p.8CI|| by means of a expansion and 
Monte Carlo simulations and in Ref. ||181|| by means of an e-expansion. The dirty boson 



system and the transition to the Bose glass phase (including the case of long-range Coulomb 



interaction) was extensively studied by Monte Carlo simulations ||182| - |185|] and Lanczos 
186| | as well as by using analytic calculations |187| . The finite- frequency 

More recently, in 



diagonalization 

properties close to the transition point were analyzed in Refs. [|188| - |190 
a series of papers, Sachdev and coworkers [|191H193|| studied non-zero temperature transport 
properties by means of a Boltzmann equation. 

A general analysis of the conductivity close to the transition can be performed based on 
scaling arguments. The frequency dependence of the conductivity criuj) has been obtained 
from its relation with the frequency dependent stiffness Ps(cc') (related to the increase of the 
free energy due to a time dependent twist) cr(a;) = Ae^ ps{—iuj) /iuj. Close to the transition 



it can be shown that the conductance obeys the scaling form ||150| , |180| , |194 | 



(37) 



where 5 measures the distance from the critical point, and where f{x,y) is a dimensionless 
scaling function. In the limit of very low temperature (compared to the frequency), the 
temperature drops out of the previous expression and the two scaling variables enter in the 
form 



(38) 



In view of the scaling behavior of the conductance, one should consider two limits uj -^T 
and uj ^ T. This point, emphasized in Refs. [|191| - |193| , is important both from a conceptual 
point of view and for a detailed comparison with experiments. The two situations correspond 
to two different experimental setups for transport measurements. In the uj <^T case, charge 
transport is governed by inelastic scattering between thermally excited carriers. In the 
opposite oj ^ T situation, collision between carriers can be neglected. In Eq.(^) the two 
limiting cases correspond to /(0,0) and /(cxo,0) respectively. It turns out that both values 
are universal but different. 

Here, we discuss the most prominent features of the conductivity close to the S-I transi- 
tion by means of the Ginzburg-Landau free energy of Eq.(^). Connections to other models 
will be given. In the linear response regime the conductivity follows from the functional 
derivatives of the partition function. In the presence of an electromagnetic potential A, the 
gradient term in the Ginzburg-Landau free energy enter in a gauge-invariant form: 



V 
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A 



By noticing that the current is the derivative of the free energy with respect to the vector 
potential and that the electric field is the time derivative of the vector potential (with a 
negative sign), the conductivity, in imaginary time, is expressed as {a{h) = x,y) 
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d rdr 



5^ \nZ 



6Aa{r,T)6A,{0) 



|a=o e 



(39) 



Using Eq. (pO[) , the longitudinal conductivity aaa = o" can be expressed in terms of two- 
and four-point Green's functions. In the absence of charge and magnetic frustration and by 



evaluating the correlators in the Gaussian approximation, one obtains [|180 

2 



(40) 



where 



G{k,ujf, 



e + + (ul 



This turns out to be the first term in a expansion [|18CI|| . The sum over the Matsubara 
frequencies can be performed by contour integration. The result is 



dz 



' dkk^QG^ik, z) \^G^{k, z) - m,G\k, z^uj) 



AttRqlu J-oo 1 - e-^^ 
where the retarded G^ and advanced G^ Green's functions are given by 

G^/^(A;, uj) = G{k, — >uj±ir]) . 
The previous expression can be evaluated in various important limits. 



(41) 



1. Zero-temperature conductivity 



Performing the /c-integral the real and imaginary parts of the conductivity are 

,2\ 
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where 
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UJr. 



UO — Ur 



UJ + UJr 



(43) 



At low frequencies the real part of the conductivity exhibits an excitation gap equal to uj^- 
In the insulating region the system, as can be deduced from the behavior of the imaginary 
part of the conductivity, behaves as an effective capacitor with 



C, 



QRqLOc 



The previous expressions can be calculated in the lowest order in and obey the scaling 
law with zu = 1. The threshold frequency vanishes at the S-I transition leading to a finite 
d.c. {lj 0) conductivity. 
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As explained in the first part of tliis section, tfiis corresponds to tfie evaluation of the scaling 
function for tu/T — > oo (the collision- free regime). Corrections to the next order in the 1/N 
expansion correct this Gaussian result by roughly 30%. 

Another powerful method for evaluating critical quantities is the e-expansion. In order 
to set up the e-expansion one should move away from two dimensions and consider systems 
with d — 1 spatial dimensions. This approach allows one to obtain the scaling form of the 
frequency dependent conductance [|181|] in d dimensions. 

In the fully frustrated case (/ = 1/2) the conductance at the S-I transition is still finite 
but with a value which is different from the / = case. It is possible to evaluate it in a 
1 /A^-expansion ||195|| and to the lowest order the conductance is twice the value of the critical 



conductance in zero field. 

a\f = 1/2) = 2a^(/ = 0) 

This factor of two is reminiscent of the superlattice structure at full frustration. There are, 
however, no fundamental reasons why this ratio should hold in general. 



2. Finite-temperature conductivity 



At low temperatures (T <^ uic) , the real part of the conductivity is given by [ |188|J191|| 
Ro 8Ro \ to 



sft ^(^) = ^Te-^-c5(^) + _^ ( 1 _ e(^2 _ n _^ 2e-^l-l/2] (45) 



The imaginary part is obtained by means of Kramers- Kronig dispersion relations 

2 



c> a(u) = _^re-^"= + uCeff (46) 



where 

^-ff = + 24^6-^-=') (47) 

The Drude peak in the real part of the conductivity arises due to a lack of dissipation 
or disorder in this model. Once electron or hole like excitations are created, they will 
propagate without damping thus leading to perfect conductivity. Although the system is a 
perfect conductor it is not a superconductor since it shows no Meissner effect. The response 
of the system to a static k-dependent magnetic field, is proportional to /c^, i.e. it vanishes 
at long wavelengths. 

The scale for the crossover to the classical behavior is set by T ~ In the high temper- 
ature limit (T ^ tUc); the real and imaginary part of the conductivity read 

«.M^.^.M + ^(i-^)e(.?-.^) (48) 
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where po ^ T and the expression for the imaginary part is vahd at frequencies much smaller 
than ojc- 

Damle and Sachdev ||191| , |192|] pointed out that since the conductance is a universal 
function of a;/T, it makes a difference which limit is taken first (either — or T — >■ 0). 
In order to study the collision-dominated regime they used a Boltzmann like approach in 
which the current is expressed in terms of distribution functions for the particle and hole- 
like excitations. By solving the appropriate Boltzmann equation (the collision term can be 
obtained by Fermi golden rule) they showed that also in the collision dominated regime the 
conductivity is a universal function at the critical point. We refer to the book by Sachdev 
for a clear and comprehensive presentation of these aspects of transport close to quantum 
critical points . 



3. Non universal behavior 

Despite the conceptual elegance of the theories predicting a universal conductance at the 
transition, the experiments on JJAs and two-dimensional superconducting films show critical 
resistivities that are different (by a factor up to ten) as compared to the predicted universal 
values. Wagenblast et al. [|1 74|] developed a theory of this non-universal behavior using 
the local-damping model. The evaluation of the dynamical conductivity proceeds along the 
same lines discussed before (see Eq. ^4|). The advanced and retarded Greens functions are 
given by 



-1 



e + k'^ — C,uj^ + ri\ujY cos — ± 2sign(c<j) sin ^ , (50) 



With increasing damping the Mott gap is smeared out. For s < 2 and low frequencies 
uj <^uJc one finds 

^ ' Rq 67re2 r(2 + 2s) ' ' ^ ^ 

The conductivity shows a power-law behavior at low frequency, where the power depends 
on the dissipation strength for s < 2. Of particular interest is the d.c. conductivity at the 
transition, which becomes a function of the strength of Ohmic damping for a > 2/3. 

This model with local damping was further explored by Dalidovich and Phillips ||196| in 



the case s = 1. In the limit of weak damping they find that dissipation leads to a leveling off 
of the d.c. conductivity at intermediate temperatures. Their estimates indicate resistance 
plateaus of the order of 10 kQ in the mK range, compatible with the experiments. These 
results seem more applicable to uniform films rather than to JJAs. In any case they offer an 
interesting explanation for the experimental observation that the critical resistance is not 
universal. 



G. One-dimensional arrays 

Josephson-junction chains have been much less investigated (both theoretically and ex- 
perimentally) as compared to two-dimensional systems and only recently the S-I transition 
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in one-dimensional samples has been measured |p,97|| . In addition to the possibility to fabri- 
cate arrays with controlled couplings, in Josephson chains the ratio of the Josephson to the 
charging energy can be varied in situ by connecting mesoscopic SQUIDS in series (as illus- 
trated in Fig.pTD. In this setup, the sample behaves as a chain of junctions with a tunable 
Josephson coupling -Ej($) = 2ii^j cos(7r$/$o) depending on the magnetic flux $ piercing 
the SQUID. By varying $, it is possible to sweep through the S-I transition while measuring 
on the same sample. 

At zero temperature and for short range Coulomb interaction the S-I transition of a 
(i-dimensional array is of the same universality class as a classical XY model in d + 1 di- 
mensions. Therefore Josephson chains should exhibit a BKT-like transition. By means of 
duality transformation it is possible to map the XY-model onto a gas of logarithmic inter- 
acting vortices ||198|| . Vortices are bound in pairs (of opposite vorticity) below the transition 
temperature and are in a plasma phase in the disordered (high temperature phase). In a 
quantum chain the relevant topological excitations, which correspond to vortices in space- 
time, are (quantum) phase slips. The mapping of a Josephson chain onto a gas of interacting 
phase slips has been performed by Bradley and Doniach ||199| ]. 

Consider, for simplicity, only the charging part of the Hamiltonian and neglect the con- 
tribution due to the junction capacitance. 



Zch = Y[ D(j)i{T)Dqi{r) exp 



Jo : , Jo 



(52) 



The summation over the winding numbers fixes the charges to be integers in units of 2e. By 
discretizing the path integral (with a time slice and performing the summation over the 
integers qi the charging contribution to the partition function can be recasted into the form 



Zch = Y[j rf0,,.^exp[-(l/8r,Eo) $]( 



(53) 



Eq.(^3]) is the Villain approximation of the XY potential ||8^ if one identifies (pi^r — 4>i,T+e 
as the dynamical variable and l/eEo as the effective coupling. The time slice e can then 
be chosen such that the coupling in space and time is isotropic (e ~ 1/ ^/SEjEq) ||200 

The XY model in space-time has a reduced coupling proportional to the ratio 



'Ej/8Eo. 

Therefore, all the known results for the classical XY model directly apply with the following 
replacement 



T 



8Eo 



(54) 



The identification of the charging energy with the effective temperature shows the anal- 
ogy between classical (thermal) and quantum (induced by charging effects) fluctuations. 
The partition function can now be expressed in terms of interacting phase slips (in the same 
fashion as in the classical where it is expressed in terms of vortices): 



Z = ^ exp 



— ^ Y P{k, uj)Go{k, uj)pi-k, -uj) 



(55) 
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where p ±1 are the "charges" associated with the occurrence of a phase shp. 
The function 

imphes that phase shps interact logarithmically in space-time. The chain undergoes a BKT 
phase transition at a critical value 



Ej/Eo ~ 2/n . 

When a Josephson coupling is larger than the critical value the phase correlator decays 
algebraically (quasi-long range order) and the chain is superconducting. Phase slips are 
bound in pairs of opposite sign and therefore they do not lead to any dissipation over a 
macroscopic region (the Josephson relation implies that the occurrence of phase slip leads 
to a voltage drop). In the opposite regime the chain is in the insulating phase. Phase slips 
are not paired and any current leads to a voltage. The correlation length is then given by 

^ ~ exp 



- {n^Ej/16Eo]y^ 

(6 ~ 1). Due to the isotropy in the space-time direction, one can now define an effective 
Coulomb gap ~ Eo^~^. As long as there is particle- hole symmetry a finite range Coulomb 
interaction does not change the universality class of the transition. A detailed analysis of 
the phase diagram, for realistic Coulomb interactions, as a function of the charge frustration 
has been performed by Odintsov [p01|| . 

The presence of dissipation (see also Section [11 ii| ) modifies the critical behavior of the 
chain. The case of a Josephson chain with Ohmic dissipation has been considered by several 
authors | [2(J2| -!2U5| by means of dual transformations and Monte Carlo simulations. The main 



conclusions of this series of works is the zero-temperature-phase diagram as a function of 
dissipation strength and Josephson coupling as shown in Fig.^. In addition to the S-I phase 
boundary there are two new phases induced entirely by dissipation: 

• for small Josephson coupling and large dissipation the chain is in a phase with local 
order. The phase difference at each junction is locked in time but the chain has no 
quasi long-range order 

• for large dissipation and large Josephson coupling there is a new type of superconduct- 
ing phase characterized by the phase slips bound in quadrupoles. 

The four different phases can be measured by considering different setups as discussed in 
Ref. mi 



We conclude this section by reviewing the experiments of Chow et al. ||197|| . The de- 



pendence of the resistance on the temperature, shown in Fig.|2^, shows a non trivial scaling 
behavior. The two set of curves (solid and dashed) refer to two chains of different length. 
While in the insulating phase the resistance increases with the number of junctions, in the 
superconducting phase the opposite trend is visible. By identifying the scale-independent 
value of Ro{T), Chow et al. were able to trace out the zero-temperature critical point 
(indicated with J* in the figure). 
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The resistance in the superconducting chains can be explained in terms of formation 
of phase shps. The flat tails in the curves are due to a finite-size effect and occur for 
temperatures of the order of the effective Coulomb gap. In this region the probability of a 
phase slip event obtained in the Coulomb gas picture presented above scales with the number 

of junctions as Nx ^ ""'^ °. Quantum phase fluctuations are suppressed by increasing the 
system size. In the insulating regime the I — V curves show Coulomb blockade with a 
threshold voltage which depends on the magnetic flux piercing the SQUID [ P06 | as shown 
in Fig.Ell Thus, there is reasonable agreement with the theory, but finite-size effects make 



a quantitative analysis difficult because of the rapidly diverging correlation lengths. 



1. ID arrays as Luttinger Liquids 

The interest in one-dimensional arrays goes further as they can be described in terms of 
the Luttinger liquid (LL) model ||207|| . The low-energy excitations of the interacting electron 



gas in one dimension are long-wavelength spin and charge oscillations, rather than fermionic 
quasi-particle excitations. Accordingly, the transport properties cannot be described in 
terms of the conventional Fermi-liquid approach. The density of states shows asymptotic 
power-law behavior at low energies. Depending on the sign of the interaction an arbitrarily 
weak barrier in a quantum wire leads to perfectly reflecting (for repulsive interactions) or 



transmitting behavior at low voltages [|208|| . It is customary to characterize this interaction 
by a parameter g such that g = 1 in the noninteracting situation while g > (<)1 in the 
attractive (repulsive) case. 



A Josephson chain seems an ideal system to explore LL correlations ||209||. In the limit 



of large Josephson coupling g = y Ej/8Eq 1; i.e., the chain behaves as an attractive LL. 
Glazman and Larkin ||210|| showed that in a certain region of parameters (close to qx = 1/2) 
between the Mott lobe and the superconducting region, there is a new intermediate phase 
which is equivalent to the chain behaving as a repulsive LL. In order to characterize this 
repulsive behavior one should consider a Josephson chain with a defect. One of the junctions 
could for example be made with a Josephson coupling much smaller than the charging 
energy |piCI| , p09|| . The different phases in the phase diagram can be characterized by the 



dependence of the Josephson current on the chain length. In the superconducting phase, 
the defect has no effect as the number of junctions increases. On the contrary, in the 
repulsive LL phase there is a strong dependence on the number of junctions. The LL phase 
can be studied also by means of Andreev tunneling spectroscopy along the lines discussed 
in Ref. |2TlJ. 



Repulsive LL behavior is also present in a Josephson ladder as discussed in Ref. ||212 
The possibility of repulsive LL behavior is related to a normal phase of interacting bosons 
at zero temperature. In one-dimensional systems such a possibility cannot be excluded and 



Monte Carlo simulations on a Josephson chain |pl3|| show a phase in which there is neither 
crystalline nor superfluid order. The existence of a normal phase has been questioned in 
Ref. [pi4|| through Density Matrix Renormalization Group of the BH model. One should note 
however that phase boundaries are non-universal and therefore the QPM and BH system 
can lead to different results. 
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H. Field-tuned transitions 



In arrays which are in the superconducting state at / = but have an Ej/Eq ratio 
close to the critical value, a magnetic field can be used to drive the array into the insulating 



state. This field-tuned transition has been considered theoretically by M.P.A. Fisher |215 



in disordered systems and has first been observed by Hebard and Palaanen pi6| , pi7|| in thin 



InOa; films. The interplay between disorder and vortex-vortex interactions plays an essential 
role. At low magnetic fields vortices at T = are pinned (by disorder) but for higher fields, 
the vortex density increases and at some critical density, vortices Bose-condense (a vortex 
superfiuid leads to an infinite resistance). By employing duality arguments this transition 
can also be thought as a Bose condensation of vortices that occurs by changing the applied 
magnetic field. 

The general characteristic of the field-tuned S-I transition is that when / is increased 
from zero, the temperature derivative of the resistance changes sign at critical values ±/c. 
Fisher's analysis |pi5|] leads to the following scaling for the resistivity tensor close to the 
field tuned transition 

where v is the exponent which controls the divergence of the correlation length at the 
transition and z is the dynamical critical exponent (with zv > 1). The resistivities are 
predicted to be universal at the transition and should satisfy the relation 

h 



PIx + Ply = ^ (57) 

These predictions were tested in Josephson arrays by the Delft pT| , |218|| and the Chalmers 
groups 
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For several values of the frustration, the resistance as a function of temperature is shown 
in Fig.|^. Below a critical value /c, the resistance decreases upon cooling down {dRo/dT > 
0). Above fc the resistance increases {dRo/dT < 0) and for low temperatures reaches a 
value that might be orders of magnitudes higher than the normal-state resistance. This sign 
change in the temperature dependence corresponds to a change in the I-V characteristics 
shown in Fig.^. For / < /c, a critical current is observed in the I-V characteristics, whereas 
above fc a charging gap develops. Note that at low temperatures, the resistance fiattens off. 
This is most likely a finite size effect involving quantum tunneling of vortices. Finite-size 
effects are expected to play a more prominent role in JJAs as compared to films because 
arrays are typically 100 cells wide. In units of the coherence length, films are much larger. 

A more detailed way of observing the field-tuned S-I transition is obtained by measuring 
the resistance versus magnetic field for different temperatures (see Fig.p7|). In the range 
< / < 1/3, the R{f) curves are very similar to the ones measured in thin films. Below the 
critical field /c=0.14 the resistance becomes smaller when the temperature is lowered and 
above fc the resistance increases. 

From the scaling analysis [pl5|] , it follows that the slopes of the R{f) curves at fc should 



follow a power-law dependence on T with power —l/{zi'). When on a double logarithmic 
plot the slopes of the R{f) curves at fc are plotted versus 1/T, one finds straight lines in 
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the temperature range 50 < T < 500 mK. /^From the reciprocal of the slope, the product 
zv can be determined. Values in JJAs range from 1.2 to 2 for the Delft data and from 
1.5 to 8.2 for the Chalmers data, in agreement with the theoretical expectations z = \ 
and u > 1. The scaling resulting from Eq.(|56D is best seen by plotting the resistance as a 
function (/ — fc)/T^^^'^ as illustrated in Fig.^ A universal function is obtained by plotting 
the resistance as a function of Ej(f — fc)/{EcT^^^^). The tails on the superconducting side 
(bottom curve) correspond to the finite-size effect mention above. 

The exponent z can also be obtained form the measurements by plotting as a function 
of the zero-field BKT transition temperature: 

The Delft-data points on the triangular arrays yield a rough estimate of z ~ 0.34 and their 
two data points on the square arrays of 2; ^ 1.4. The Chalmers data provide a more accurate 
fit yielding z = 1.05, in good agreement with the theoretical expectation. 

Measurements on different thin films show that the resistance right at the transition is of 
the order of Rq but measurements are not conclusive regarding the universality. In arrays, 
this resistance is again of order Rq, yet in different arrays it varies between 1.6 and 12.5 kQ. 
The Chalmers group has also measured the Hall resistance in order to check the validity of 
Eq .|57|. For two ar rays, the Hall resistance at the critical point is of the order of 30 Q, but 
again ^pl^^ + pl y is not a universal quantity. 

A new feature introduced by JJAs is the existence of field-tuned transitions near com- 
mensurate values of the applied field, i.e., at fcomm ± fc ||53| , |218|| . Studying the R{f) curves 



of JJAs in more detail, critical behavior is not only seen around / = 0, but also around 
/ = ±1/4, ±1/3, ±1/2, ±2/3, and ±3/4. For each commensurate /-value zu can be de- 
termined as described above and the values of zu are close to one. The sample-dependent 
critical resistances are of the order of a few kfl. Calculations on the Bose-Hubbard model in 



a magnetic field [ 114 | show that the product zu ai f = 1/2 is close to 1 in agreement with 
the measurement. 



III. QUANTUM VORTEX DYNAMICS 

A vortex (antivortex) is a topological excitation. When going around its center in a closed 
loop, the phases of the islands wind up to 27r (-27r). Vortices have extensively been studied 
in classical arrays where they both determine both the phase diagram and the dynamical 
properties. To a large extent, this still holds for quantum arrays where the interplay between 
vortex and charge dynamics plays the central role. The field tuned transition for instance 
can be understood as a Bose condensation of vortices or charges and the S-I transition in 
zero field can be analyzed using the duality between charge and vortex excitations. In this 
section we show that vortices characterize the quantum dynamics of arrays as well. We 
concentrate on the superconducting side of the S-I transition where vortices are well defined 
excitations. 

As discussed before, quantum arrays have flux penetration depths that are larger than 
the array sizes. The magnetic field is almost uniform over the whole array area indicating 
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that there is not one flux quantum in particular cell. f\ 

The essential aspect of vortices in junction arrays is therefore not the flux, but the 
distribution of phases. The phase conflguration of a vortex (shown in Fig.|29|) in a large 2D 
array can be approximated by the following analytical expression 



± arctan ( — — — ) (58) 



where the site i has coordinates = and the vortex center is placed at (r = [x, y]). 

The ± sign refers to the vortex (antivortex) conflguration. For most purposes Eq.(|58D is 
accurate even very close to the vortex center. As we will see, the approximate arc-tan 
solution is very useful since it allows to express the action in terms of the coordinates of 
the vortex center r only (instead of in term of all the phases). This appears to be a reliable 
description as long as the vortex can be considered as a rigid body. In most of the cases we 
discuss in this review, this turns out to be a good approximation. 

An important property of vortices is that supercurrents around them fall off inversely 
proportional to the distance r from their core. Vortex-vortex interactions therefore have 
a long-distance character as they are proportional to Inr. Eq . (|58|) is the solution for a 
single vortex in an inflnite system. In flnite systems, vortices interact with boundaries. The 
interaction of vortices with the open edges can be viewed as the attraction of a vortex with 
an image antivortex outside the array. Superconducting banks repel vortices; the interaction 
with these edges can be viewed as the one with an image vortex (of the same sign) outside 
the array. The interaction between boundaries indicates that especially in small arrays the 
approximation given in Eq.(^) is no longer valid. Numerical calculations are then used to 
extract the quasi-static phase conflguration around a vortex. 

Experimentally, single- vortex dynamics is studied by applying a small magnetic fleld (low 
vortex densities) and performing transport measurements. On the theoretical side, both 
numerical simulations and phenomenological models which lump the collective dynamics of 
the phases into the description of the motion of the vortex center, have been investigated. 
In this chapter, we flrst derive the classical equation of motion and show that vortices 
in underdamped arrays can be viewed as massive, point-like particles. We then continue 
with the quantum corrections to the equation of motion and discuss their consequences. 
First, there is a renormalization of the vortex parameters (like the mass, damping, ...) 
due to quantum fluctuations. In addition there is a class of new phenomena which arise 
from the quantum dynamics of vortices which are treated in Sec. [1IID| and |III E| . They 



include macroscopic quantum tunneling of vortices in 2D arrays, quantum interference in 
a hexagon-shaped array, Bloch oscillations of vortices in the periodic lattice potential and 



^Self-field effects may play a role in classical arrays since critical currents are substantially larger. 
Generally, speaking self-field effects manifest themselves in two ways [21£]. First, there are self- 
inductance effects which are short-ranged and caused by the self-inductance of the cell loop. It 
turns out that the cell-to-cell energy barrier is dominated by these short-range interactions. Second, 
there are mutual-inductive effects which have a longer range. For example, the current distribution 
around a vortex changes from an exponential fall-off for self-inductances to an algebraic fall-off when 
the mutual inductances between all cell pairs are included. 
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vortex localization in quasi one- dimensional arrays. In some experiments, evidence has been 
found that the (independent) single-vortex picture breaks down. In these cases, we briefly 
comment on the influence of vortex-vortex interactions. 



A. Classical equation of motion 

In this section, we review the steps to derive the equation of motion for a vortex. We 
analyze all contributions, i.e. its inertia, its dissipation, the external potential and the 
applied forces. To keep the notation simple we suppose, for the moment that the vortex 
moves along a given (say x) direction with an average vortex velocity Vx- 

Vortex mass - Moving vortices lead to phase changes across junctions and they therefore 
contribute to an electric energy. In principle all capacitances contribute, but as discussed 
before the main effect comes from the junction capacitance so that Ech = l/2Z]<ij> CVij)- 
(The contribution due to Cq will be discussed in Section [ill C] .) In a quasi-static approach 
this sum is calculated by comparing the phase differences across each junction at times t 
and t + 1/vx (in units of the lattice constant): A0jj = (pijit + 1/vx) — 4>ij(t)- The electric 
energy then acts like a kinetic energy term and the proportionality factor deflnes the vortex 



mass 



= ^M.vl (59) 



with 



Mv = ^ E(A</>..)^ (60) 

<ij> 



8Ec 



The problem of calculating the vortex mass is now reduced to flnding the phase differences 
across junctions at times t and t + l/vx- Note, that Eq.(^D]) can be applied to various array 
geometries if the phase conflguration around a vortex is known. 



In large 2D Josephson arrays the arctan form given in Eq. (|58|) [|224|| is used to evaluate 
the sum in Eq. (|60|). With the assumption that this arctan-phase conflguration remains the 
same when the vortex moves through the array, numerical evaluation of the phase differences 
in a large 2D square array yields the Eckern-Schmid value of the vortex mass 

M, = ^Ec' . (61) 

In this calculation roughly half of the vortex mass is due to the junction the vortex crosses; 
the other half comes from all the other junctions in the array. For a triangular 2D array, 
a similar calculation can be done and the vortex mass is twice the mass of a square array. 
For typical arrays with C = 1 fF and = 10 /xm^, the vortex mass is 500 times smaller 
than the electron mass. This small value already indicates that quantum effects are likely 
to occur. 

It has been shown ||225|| that near array edges the vortex mass vanishes when it approaches 
a free boundary of the array. These boundary effects are, however, negligible if the vortex is 
a few lattice spacings away from the edge. One can also include self-fleld effects. Currents 
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now extend over a distance of the penetration depth A_l from the vortex center so that the 
arctan approximation can no longer be used. As the vortex is effectively reduced in size, the 
sum of the Vj's can be restricted to those junctions which are A_l from the vortex center. The 



result is a smaller vortex mass and its decrease with decreasing A_l is given in Ref . ||226|| . The 
vortex mass is dependent on the Ej/ Ec ratio as well. Quantum corrections to the mass, on 
approaching the S-I transition ||227|| are discussed in Section [III C| . To a good approximation 



it has the value given in Eq. ( pT] ) for arrays that are not in the critical region Ej ~ Ec- 

The vortex mass can also be calculated in geometries other than two dimensional arrays. 
In a purely ID array (A^ junctions in parallel connected by two superconducting leads), the 
vortex phase configuration is given by 

(pi = 4 arctan[exp((xj — x)/Aj)], 

where x denotes the position of the vortex center. For Aj < A^, the vortex has a kink- 
like shape, which extends over a distance of the order of Aj; for Aj > N the vortex is 
spread out equally over the whole system with (px+i — 4>x ~ 2tt/N. In this latter regime, 
A(f)t+i/v,{x + l,x) - A0^+i,^(t) = 27i/N and = h'^/{SEcN^). For 1 < Aj < N the sum 
can be computed numerically or in a continuum approximation. The sum over the phase 



differences squared is equal to 8/Aj ||228|| and hence 



Mj^ = ^ Ec^ for 1< Aj < A^ . (62) 
Aj 

In other geometries, the phase configuration of a vortex is not exactly known. It can in 
principle be calculated. However, since a substantial contribution comes from the junction 
that the vortex crosses, Eq. (^) can be used as an estimate for the quasi-static vortex mass 
in these cases. 

The notion that the concept of the vortex mass is not new in the theory of superfluids. 
It has been discussed extensively for type-II superconductors and superfluid He ||229| - p31| . 



Dissipation - As the vortex is a macroscopic object, it couples to the environment and 
experiences dissipation. Quantum arrays generally have junctions that are underdamped, 
i.e., have a McCumber parameter /?c > 1- One might therefore expect that vortex motion is 
underdamped as well. 

In the simplest approximation one can assume that a moving vortex experience a viscous 
drag force characterized by a viscous coefficient rj. In a Bardeen-Stephen like model t] is 
calculated using the following argument. The total power loss is the sum of all the resistive 
losses in the junctions. Assuming to be identical for all junctions, the sum in the total 
power is the same as in the calculation of the vortex mass. For example. 



2R, 



for a large square 2D array. Here, Rg is the effective voltage-bias resistance, i.e., the effective 
shunt resistance of each junction. At low temperatures Re is the subgap resistance which is 
many orders of magnitude larger than the normal state resistance, indicating the vortices can 
move through the medium with negligible damping. However, the simple model presented 
here does not take into account other sources of dissipation like the coupling to the low lying 
modes of the array (spin- waves) or quasi-particle tunneling (see next sections). 
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Lattice potential - The total Josephson energy associated to a vortex configuration (cal- 
culated for example by means of Eq.(^)) depends on the vortex position. The energy has 
a minimum value when the vortex is in the middle of a cell; the maximum value is reached 
when the vortex right on top of a junction. Neglecting vortex-vortex interaction and the 
infiuence of the array edges, vortices are only subject to a periodic lattice potential: 

Uy{x) = ^'yEjsm{27!-x) . 

Here, 7 is the energy barrier in units of Ej a vortex has to overcome when moving from 
one cell to the next. In large 2D arrays with no self fields, 7 = 0.2 in a square geometry. 



In a triangular arrays the barrier is about a factor five lower, 7 = 0.043 [|224|| . Inclusion 



of self-field effects can be done and 7 increases dramatically for A_l < 1 ||219|| . In contrast, 
there is no energy barrier in ID arrays if Aj > 1 |p32|] . 



Equation of motion for a single vortex - A vortex in a Josephson array moves under in- 
fiuence of a Lorentz force $0-^ in a direction perpendicular to the current fiow. The phe- 
nomenological damping term and the periodic lattice potential U{x) provide additional 
forces. 

Gathering all the ingredients discussed so the equation of motion can be written as: 

M„r + r]r = -W^U^ - <l>oz x I (63) 
where r is the vortex position, I is the applied current per junction and z is the unit vector 



perpendicular to the array |233 



The dynamics can be visualized as that of a massive particle moving in a washboard 
potential analogous to the dynamics of a single junction in the RCS J model. For the junction 
problem, the motion is in artificial 0-space; for the vortex problem the motion is in real space. 
The mapping is exact for 27rx — ^ indicating that vortices in arrays produce the same 
dynamics as a single junction with a critical current of 7/c/2 per junction, a McCumber 
parameter /3c,t, = 7/3c and a plasma frequency ujp^^ = y/^ujp 



.p. 



Numerical evidence for massive vortices has been found by Hagenaars et al. ||234|| . Their 
data show that a vortex may be refiected at an array edge, thereby changing its sign (i.e., 
it becomes an antivortex). This behavior can be qualitatively understood within the model 
of a massive vortex interacting logarithmically with the image vortices outside the array. In 
the same paper the authors also note that the way in which vortex inertia manifests itself 
depends on the dynamical situation considered. 

In the next three subsections, we first summarize the experimental details on classical 
arrays with (3c^v > 1 and then discuss two phenomena that are not included in Eq. (|63D. First, 
experiments show that vortices in highly underdamped arrays experience more damping 
than can be expected from the simple approach we followed above. As it turns out coupling 
to spin waves becomes dominant in these arrays. Second, we review the theoretical and 
experimental results on the Hall effect. 



1. Experiments on classical, underdamped arrays 



In Fig.^ a typical example ||235| , P36|| of a current-voltage I-V characteristic of a classical 
2D array is shown. The applied magnetic field corresponds to / = 0.1. At low temperatures. 
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hysteresis near the depinning current indicates that [3c,v > 1, consistent with the existence 
] of a mass term in the equation of motion. The depinning current itself is close to the 
expected value of (7/2)/^ = 0.1/c per junction. From the analogy with the single junction 
problem, a RCSJ-like I-V characteristic is expected. This is generally not observed. The 
I-V curves, instead, show a slight bending in the direction of the voltage axis opposite to 
what is expected from the single-junction analogy. This is also seen in simulations on the 
properties of a single vortex in a 2D array with periodic boundaries | ]237| . Their numerical 
data points at a nonlinear viscous damping of the form: 

A 

V 



where the A and B depend only on the McCumber parameter of the junctions. 

For currents well above depinning (above 50 /xA and not visible in Fig.|3y), the flux- flow 
state becomes unstable. The I-V enters a row-switched state [p.05|| where rows of junctions 



across the whole array width start to oscillate coherently |p38|] : all phases rotate continuously 
in time with a phase shift between them. In this regime, a description of the array in terms 
of vortex motion is no longer appropriate. 

One should realize that the I-V curve was recorded at an applied magnetic flux of 0.1 $0 
per cell (/ = 0.1). Thus, on average there is approximately one vortex per 1/f cells so 
that the distance between vortices is only three cells. At such a short distance, vortices will 
interact with each other. The influence of these vortex-vortex interactions on the measured 
I-V characteristics is not known in detail. 



2. Spin-wave damping 



As shown in Fig.^ experimental I-V characteristics in the flux-flow regime are generally 
more or less straight lines. Neglecting the influence of the pinning potential U^, Eq.(|53D 
indicates that for high bias the slope of this line should approach a conductance value 
corresponding to l/(2/i?e) per junction. The I-V curves therefore provide a way to estimate 
i?e in the regime where vortices are driven with relatively large currents. 

A systematic study on highly underdamped arrays has been performed by the Delft 
group [p39|| . The surprising result is that for the most underdamped arrays Re is much 
lower than the normal-state resistance. Such a low resistance cannot be explained by the 
Bardeen- Stephen model. Apparently, vortices, when driven with a large current, experience 
more damping than can be explained by ohmic dissipation alone. A similar conclusion was 
drawn by Tighe et al. [ 240 |, who concluded that in their underdamped arrays vortices moved 
in an overdamped manner. 

Several authors have suggested the possibility that energy can be lost in the wake of the 
moving vortex [ p35|j239|J240|| . The effective viscosity due to coupling to spin- waves can be 



calculated in a semi-quantitative model ||239|| using the following argument: The oscillating 
part of the junction is modeled by an Lj — C circuit and the voltage drop across it is 
V = $o'f^^0/(2vr) which occurs for a time interval of the order of v~^. As a response to 
this voltage step the phase difference starts to oscillate (note that at the same time the 
average phase advances when the vortex moves across the junction). Equating the total 
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power dissipated in the Lj — C circuit to rjswv'^, the result for the effective viscosity due to 
the plasma oscillations (in square two-dimensional arrays) is 

1 $2 1 



When comparing this viscosity coefficient to the Bardeen-Stephen viscosity coefficient, 



vv ,/lJ/c 

one sees that the more underdamped the arrays are, the more dominant the damping due 
to energy lost in the wake of the vortex becomes. 

These observations have been confirmed by more systematic calculations discussed in 
Section [III Q . In particular it is possible to obtain a self-consistent picture of vortex dynamics 



which includes the interaction with its environment. There are two main advantages for 
discussing the dynamics from this perspective. Firstly it is possible to evaluate quantum 
corrections to the classical equation of motion and secondly it is possible to analyze the 
quantum dynamics in detail. 



3. The Hall Effect 

In addition to the Lorentz force which is due to the external current, a vortex is subject 
to a Magnus force which is transverse to vortex velocity. The study of the Magnus force in 
superfiuids has a very long history. A detailed discussion is outside the scope of this review 
(see Ref. ||241|| and references therein). In Josephson arrays, in presence of a gate to the 



ground plane, particle-hole symmetry is broken. A vortex feels a Magnus force [p42|j243 
given by 

F = Q,$oz X r . (64) 

Here we assumed for simplicity a homogeneous gate charge. As a result of the combined 
effect of the Magnus force and the Lorentz force, the vortices move at a certain angle, the 
Hall angle, with respect to the current. Its measurement yields informations on the different 
dissipation sources in the system. Combining all the terms the equation of motion in the 
stationary limit the vortex moves at a constant velocity v = [fa;,^?;] = [f cos 6'^, f sin 6*^;^] 
obeying the following equations 

rjv cos 6h = /y$o — QxV^o sin 9^ 

rjv sin = Q^v^q cos 9h (65) 



lead to the resistance tensor 

Rxx Rx 



1 + {Qx%hY 



p _ p _ QMh^ /RRN 
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and to the Hall angle 6h 



tan 6 



H 



(67) 



The main consequence of the previous results is that the Hall effect should be larger in low 
resistance samples. 

In experiments on classical Josephson arrays the Hall angle is usually found to be very 
small (see e.g. Ref. p44| ). Samples are usually characterized by random offset charges and 
as a result the Magnus force averages to approximately zero. However, up to now there is no 
general agreement on this explanation. From a theoretical point of view there are questions 
related to the derivation of Eq.(0) from first principles. 

In Ref. |p42 the Magnus force was obtained from the QPM, implying that only the 
external charge enters in determining the Hall angle. A reexamination of the problem by 



Makhlin and Volovik ||245|| related the (apparent) absence of the Hall angle to the near exact 
cancellation of the Magnus force with the spectral-fiow force. On deriving the effective 
action from the BCS Hamiltonian Volovik ||246|| shows that the offset charges, contributing 
to the Hall angle, have two different physical origins. In addition to the one stemming from 
the coupling to the ground plane, there is an additional contribution which depends on the 
particle-hole asymmetry of the spectrum. This latter term is of the order of the small factor 
{A/EfY, with Ep being the Fermi energy. This confirms the expectation that the Hall 



angle should be small in Josephson arrays [|247 



In quantum Josephson arrays Hall measurements have been performed by the Chalmers 
group ||161| , p48|| . Their results are shown in Fig.|3ll. The transverse resistance is odd in 
the magnetic field. Combining the results for the longitudinal and transverse part the Hall 
angle can be extracted. Comparing the results with the classical results [p44|| , the Chalmers 



experiments indicate a larger Hall angle. The only apparent difference is the smaller ratio 
Ej/Ec and this is consistent with the theoretical expectation that the offset charges are 
responsible for the Hall effect. It is reasonable to expect that these have a negligible effect 
on approaching the classical limit. 

Finally, it is worth mentioning that the field-tuned transition discussed previously also 
manifests itself in R^y. An interesting feature which still remains unexplained is that Rxx 
andRxy are related by the following empirical relation 



Rxy ~ 



dRxc 



(68) 



similarly to what happens in the Quantum Hall effect. 



B. Ballistic vortex motion 

Besides the experimental verification of the mass term in the equation of motion [ ^35|j24CI|| 
a considerable interest was focused on the direct observation of the ballistic motion. Ballistic 



vortex motion has not only been observed in long continuous junctions [p49|] , where energy 
barriers for cell-to-cell motion and spin- wave coupling are absent, but also in discrete ID 
arrays ||25CI| - p52|] and in 2D aluminum arrays ||253|| . The idea goes as follows: If vortices are 
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massive particles, they should keep on moving if the current is turned off. In an experiment, 
this concept can be realized by accelerating vortices up to a high velocity Vq so that their 
kinetic energy is much larger than the lattice potential. With Eq.(|63D one finds that Vq ~ 
^ol/t] if one neglects the lattice potential. Then, these fast-moving vortices can be launched 
into a force-free environment where voltages probes can be used to detect their path through 
this region. 

The observation of ballistic motion is only possible in a velocity window {vmin < v < 
Vmax) bound from below by the presence of the pinning potential and from above by the 
various damping mechanisms which set in at high velocities. The criterion E^in > Ep^t 
translates into a lower bound for the vortex velocity required to observe ballistic motion 

Vmin = • (69) 

71 

Note, that for a ID system with Aj > 1, 7 so that the minimum vortex velocity is 
small. 

The vortex velocity cannot be chosen arbitrarily large. Fast moving vortices can trigger 
row switching in the array ||105| , |254|| . Simulations ||255|| indicate that in 2D arrays the vortex 



velocity must be limited to < Up. Another limitation comes from coupling to spin-waves. 
In 2D arrays, there is a threshold vortex velocity below which this coupling is weak. It has 



been shown ||256|J257|| that a moving vortex only couples to spin- waves above v^ax ~ O.luj. 



V 



The requirement that 



77' 



uo. 



p 



TC 

indicates that ballistic motion is possible in triangular arrays just above depinning. (Quan- 
tum fluctuations make the window larger; details are discussed in the next section.) 

Let us analyze the ballistic motion using Eq. (|63D . With no current applied and for 
vortices launched at high enough velocity (to ignore pinning effects), the equation of motion 
reduces to 

Myi)^ + rjVx = , 

indicating that the vortex velocity decreases exponentially in time as Vx = f (0) exp[— M^yt/r^]. 
A mean free vortex path can then be defined as 

v{0)M, I 

^free = = TT — Pc , (70) 

for a square 2D array. The factor I/ttIc is typically of order 0.1 so that at high temperatures 
with Re = Rn, ^free ~ 1- At low tcmpcratures with R^ ^ R^ (the corresponding Pc can 
be high as 10^), Xjree ^ 1- 



The experiment with 2D arrays was performed by van der Zant et al. using the 
sample configuration shown in Fig.^. It consists of two 2D arrays which are connected by 
a narrow channel of 20 cells long and 7 cells wide. Superconducting banks on both side of 
the channel confine the vortices in the channel and, in order to reduce the influence of the 
lattice potential, the arrays and the channel were made in a triangular geometry. 

In the array on the left hand side, vortices generated by a small magnetic field, are 
accelerated up to a high velocity. Some of these high-energetic vortices will enter the channel 
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and will then be launched into the detector array (array on the right hand side). There is 
no driving current applied to this array. A set of voltage probes around the detector array 
is used to detect the places where vortices leave the force- free environment. (The voltage 
measured across two probes is proportional to the number of vortices passing the probes per 
unit time.) The results of the local voltage drops as a function of temperature it is shown in 



Fig.^. At high temperatures, vortices move diffusively and voltages are observed between 
all voltage probes. As a consequence, the voltage across the probes V3 and VIO is much 
smaller than the channel voltage. 

At low temperatures, the voltage measured between the two probes situated just opposite 
to the channel is almost equal to the channel voltage. Vortices cross the second array in 
a narrow beam (see Fig.^2]). This ballistic vortex motion is observed for T < 500 mK, for 
small applied magnetic fields (0.01 < / < 0.025) and for currents just above depinning. 
For high magnetic fields, vortex-vortex interactions start to play a role when more than one 
vortex is in the channel at the same time and for too high currents coupling to spin-waves 
probably starts to play a role. 



C. Effective single vortex action 

In this Section we derive a more general approach to describe the single vortex dynamics 
which incorporates the coupling to spin-waves and which is also valid in the quantum regime 
where Ej ^ Ec | |13U| , ^27|j258|| . After introducing the effective vortex action, we first consider 
the classical limit and show that all known results can be recovered. We then proceed to 
the quantum regime. The vortex mass is calculated as well as the velocity above which 
spin-wave damping starts to become effective. 

As discussed in Appendix the effective action for a single vortex is given by 

Seff = 1 J dTdT'r\r)NU[v{r) - r(r'), r - r']f^(r'), 

a,b=x,y 

Mat = E Vae(r(r) - r,){q,{r)qk{r'))VM^k - r(r')) (71) 

jk 

(0(r) = arctan{y / x)) . Thus, vortex dynamics is governed by the charge-charge correlation, 
which depends on the full coupled charge- vortex gas. The effective action Eq.(|7TD describes 
dynamical vortex properties in the whole superconducting region and is therefore a good 
starting point for the investigation of vortex properties down to the S-I transition. The 
cumulant expansion that leads to Eq. (|7TD is correct in the Ej 3> Ec limit where the charges 
can be considered as continuous variables and where the vortex fluctuations can be disre- 
garded. In general the average defined in Eq. (|El| ) is far from being gaussian so that one may 
argue that higher order cumulants should be considered. Nevertheless nothing prevents us 
to analyze Eq.(^) keeping in mind that a full description of vortex motion may require the 
analysis of a dynamical equation that contains also terms proportional to higher powers of 
the vortex velocity. 

The expression given in Eq. ( |7TD reproduce the known results in the classical limit where 
Ej ^ Ec- In this region of the phase diagram the charges may be considered to be continu- 
ous variables and vortex fluctuations may be neglected so that the charge-charge correlation 
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reads 



and 



Co 1 + A2A;2 • 

The spin-wave dispersion is described by ujk- It is optical, i.e. 



for long 



range Coulomb interactions, whereas for on-site interactions we have = Upk. Here 



p — y Ae^Ej/Co is the plasma frequency for the case of on-site Coulomb interactions. 
The action ( |7TD reduces to that of a free particle in the limit of small velocities r(r). The 
corresponding adiabatic vortex mass My is 



which reduces in the classical limit to 



/3 



TT 



4E, 



rfrM,.,(0,r) 



+ ^ln(L) 



c 



4e2 



with L the array dimension. Thus both Cq and C yield a contribution to the mass. The 
self-capacitance contribution depends on the system size L. For generic sample sizes and 
capacitance ratio's the size-dependent contribution is smaller than the Eckern-Schmid mass 
(the value in Eq. (0)). The effect of a uniform background charge on the vortex mass was 
considered by Luciano et al. |P5S ,260|. The frustration of charging leads to a renormalization 
of the mass towards the classical value. 

The spin-wave damping that a moving vortex experiences may also be calculated from 
Eq. ([7TD . Varying the vortex coordinate r"(r) in Eq.([7T]) yields the equation of motion 



27Teabh/Ic = dr J cLt' Mab{r{r) 



(72) 



~yy 



0,e 



xy 



— e 



yx 



1). Its constant velocity solutions in the presence of an external 
current determine the non=linear relation between driving current and vortex velocity (i.e. 
the current- voltage characteristics), once the charge-charge correlation is analytically con- 
tinued to real frequencies (i.e., if iui^ u + |P56|| . The relevant information is contained 
in the real part of Eq.(|7^), which reads in Fourier components and for a constant vortex 
velocity 



(73) 



The delta functions express the spin-wave dispersion (from the analytic continuation of 
the charge-charge correlation) and the vortex dispersion respectively. The overlap integral 
determines the amount of dissipation a moving vortex suffers from coupling to spin-waves. 
Adopting the smooth momentum integration cut-off, introduced in Ref. [ p21| 
in the classical limit the results of Refs. |p56| , p57| . 



one recovers 
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While the static damping is zero for vortex velocities below threshold (which implies the 
possibility of ballistic motion), a dynamical friction due to the coupling to the plasma oscil- 
lations is always present for frequencies higher than a given frequency threshold |p61|] . The 
latter contribution approaches to zero when the velocity increases to the threshold velocity. 
However, radiative dissipation of the vortex affects the threshold for ballistic motion. What 
is important in this analysis is that by changing the frequency of the applied current, one 
is able to extract the domain of validity where a vortex can be defined as a macroscopic 
object. 

Inclusion of quantum effects contributes to the opening of a more robust velocity inter- 
val where ballistic motion can be observed. When the ratio Ej/Eq decreases the charge- 
charge correlation must be calculated beyond the classical approximation: The discreteness 
of charge transfer has to be taken into account, which is in particular important at short 
distances. For long range Coulomb interactions and in the absence of vortex fluctuations 
the charge-charge correlation function may be rewritten as 

{qq)k,., = ^^%2^ ojl = u;l + A'K^Ejek\ (74) 



, i2 _|_ ,-,2 ' 



where the correlation length is 




vre 



and the constant c is of order one. The S-I phase transition takes place at Ej/Ec = l/vr^. 
Thus, without vortex fluctuations the phase transition is at a smaller Ej/Eq value than 
the 2/7r^ that follows from a duality argument ||9l|. The spin- wave dispersion is affected at 



small distances (large k) and the mass is now given by: 



TT 



M = -In 

IQEce 



1 + 



In the limit of small ^ the Eckern-Schmid mass is recovered. An extrapolation to the S-I 
transition where —>■ oo yields a mass that vanishes at the transition. 

With the charge charge correlation given in Eq.(|7^ we may calculate the spin- wave 
damping of vortex motion due to the coupling to spin- waves beyond the classical limit. 
Replacing Uk by ujk in Eq. (|75D , the overlap integral over the delta functions only contributes 
for vortex velocities that are higher than a threshold velocity 



~ ^\jSEjEc 

By taking into account quantum effects the spin- wave spectrum enlarges the velocity window 
in which vortices move over the lattice potential without emitting spin-wave. An extrapola- 
tion to the S-I transition yields a diverging threshold Velocity so that vortices and spin-waves 
decouple. Note, that the outcome of this calculation also has consequences for the classi- 
cal equation of motion. For instance, it shows that coupling to spin-wave dissipation is 
reduced for velocities < O.ScUp. This value is a factor five higher than expected from classical 



considerations 256,257 
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Spin-wave damping is not the only source of dissipation. Vortices, in their motion, can 
excite quasi-particles as well if the local voltage drop (due to the finite velocity) exceeds the 
quasi-particle gap. The effect of quasi-particles damping on vortex motion was considered 
in Ref. ||262|| . In this case the equation of motion takes the form 



r]r - r/A dt' sinc[2A(t - t')] 



7T\r{t) - r(t') 



-<l>oz X I 



(75) 



where sincx = (2/7rx)sinx. Despite the nonlinear form of the damping kernel, Choi et 
al. p62|| showed that the frictional force on a vortex is linear in the vortex velocity for 
any practical purpose ( in particular in the long-time and long- wavelength limit, where the 
semiclassical equation of motion is mostly concerned). 

D. Quantum vortices 



If vortices are massive particles that move ballistically, one can think of them as quantum 
mechanical objects. Like an electron, a vortex in a periodic potential will have a Bloch wave 
function with momentum p = hk and thus a wavelength of h/vM^. At present, many experi- 
ments have verified the concept of a quantum vortex ||5T| , |53| , |235| , |240| , |263| - |265|| . In this chapter, 
we discuss three examples: macroscopic quantum tunneling of vortices [ ^|53| , p35ip40|] , the 



observation of vortex interference in a hexagon- shaped array [p66|| and Bloch oscillations of 
vortices in the periodic lattice potential ||267 . 



1. Macroscopic quantum tunneling of vortices 



In a classical description vortices oscillate in the minima of the washboard potential 
with frequency ujp^v In quantum arrays, these oscillations are quantized. To estimate when 
quantum fluctuations in the vortex position become important, we compare the zero-point 
energy ^hup^ = ^^/8^yEJEc to the energy barrier f/bar = jEj. The two energies are equal 

if Ej/Eq = ^2/7. In this quantum vortex regime, the zero-point fluctuations are large 
enough to allow for quantum tunneling of vortices. 



In Fig.^ the resistance per junction (linear response) as a function of temperature is 
given for two square arrays |]^; in the classical (a) and in the quantum regime (b). The 
resistance of the classical array decreases exponentially all the way down to the lowest 
temperatures. The slopes define the barrier for this thermally activated process: R oc 
exp [—Uha.r/kBT] = exp [—PjEj]. In contrast, the resistance of the quantum array levels off 
at a temperature Tcr below which it remains constant. We denote this constant value with 
Re- Above Tcr, again thermally activated behavior is observed. Similar resistance curves 



have been reported by the Chalmers group [Q. 

A first estimate of the tunnel rates and of Rc can obtained from the analogy with the 



single-junction problem. In the moderate damping regime [268 



Rc ^ UORQfVSe- 



(76) 



where the action S is given by 
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s 



/0.95 



Ec 



1 + 



0.87 



(77) 



The s^EjjEc dependence of the critical resistance, imphed by the previous equation, has 
indeed been reported by the Chalmers group [p3 |. 

An estimate for T^r can be obtained by equating S to ^EjjTcr- Neglecting the term with 



/?ci,, the result is: 



(78) 



With Ej of the order of Ec-, the inverse normalized critical temperature (Ej/T^r) is typically 
somewhat larger than 2.5 in agreement with the data. 

Comparing the measured values of Rc with the estimates given in Eq.([7^) and Eq. ([T?]) , 
the tunnel rates in the measurements are lower than expected even when taking Rj^j as the 
resistance determining f]c,v A smaller R^ is consistent with a single vortex model in which 
the vortex mass is an order of magnitude larger than the one calculated in the quasi-static 
approximation. It is likely that vortices do not move as rigid objects and calculations have 



shown that the dynamic band mass of a vortex can be an order of magnitude larger |269 



However, considering this uncertainty, no definite conclusion about the validity of the single- 
junction model can be drawn for the observed flattening of the resistance. Other models 
like collective tunneling cannot be excluded. 

A surprising result is that the array (a) in Fig.p3 does not show any signature of quantum 
tunneling. Our simple argument given above indicates that for this array the zero-point 
energy is of order t/bar- The absence of quantum tunneling is explained by the fact that 
in Fig.B4| the measured energy barriers are of order Ej, instead of 0.2Ej. The Delft group 



has reported [^] a systematic increase of the measured energy barrier in the range 2 < 
Ej/Ec < 20. This increase is not yet understood. 

The theoretical analysis of inter-site vortex tunneling in Josephson arrays was first formu- 
lated by Korshunov [|270| , p71|| who evaluated the instanton action Sinst associated to vortex 
tunneling between adjacent plaquettes. Smst, related to a hop from one plaquette to a 
neighboring one, determines tunnel rates and the depinning current. It can be obtained in 
the language of the Coupled Coulomb gas approach by evaluating the action associated 
with the trajectory 

for a hop from x,t ^ x + 1,t + r^, with the result (see Eq.([7lD) 



inst 



:M,,(0,0). 



In the limit of large Josephson coupling one recovers all known results, i.e. for general 
capacitance matrix 



S.. 



inst 



ttEj 
4uj„ 




(79) 
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It reduces to Sinst = T^^^'^Ej/Aup and Sinst = '^'^Ej/2ujp for C = and Cq = respectively. 

Korshunov pointed out that instanton-instanton interaction cannot be neglected. Vortex 
tunneling is incoherent in the temperature range Ej ^ T ^ y/EjEo. In this case the 
tunneling probability W is given by 

W ~ ^ e2nln2T/Ej 

2^3/2 

up to the crossover temperature T ~ y/EjE^ where the activated behavior takes place. 
Dissipation associated to vortex tunneling was discussed by loffe and Narozhny ||272| . 



Since the the time associated to vortex tunneling is slow compared to A~ , the dissipation 
which accompanies this process arises from rare processes when a vortex excites a quasi- 
particle above the gap. These authors find that this source of dissipation can be significant 
even in the adiabatic limit. 



2. Vortex interference: the Aharonov-Casher effect 



In 1984 Aharonov and Casher |p73 ] studied the interference of particles with a magnetic 



moment moving around a line charge. This AC effect is the dual of the Aharonov-Bohm 
effect p74|| that describes quantum interference of charged particles moving around magnetic 



flux. The AC effect has first been observed using neutron beams in Ref. ||275|| . The concept 



of vortex-interference in superconductors has been introduced by Reznik and Aharonov ||276| 



and their ideas have been adapted to ring-shaped Josephson arrays by van Wees ||277|| and 



by Orlando and Delin ||278|| . Although there are similarities with the conventional AC effect, 
there are important differences. In JJAs vortices do not carry a flux unlike the Abrikosov 
vortices in a bulk superconductor. Moreover, as already stressed, JJAs form an artificial 2D 



space. The observation of the AC effect for vortices has been reported by Elion et al. ||266 



only two years after the first theory papers ||277|,|278|| appeared 



In the array of Elion et al, vortices follow trajectories indicated as dotted lines in the 
Fig.p^. The sample consists of a hexagon- shaped array with six triangular cells. Large- 
area junctions couple the hexagon to superconducting banks so that only two paths for 
vortex motion are possible. The large- area junctions confine vortices to the hexagon, but 
the coupling to the superconducting banks is not so strong that phases of the islands are 
set by the banks. A gate controls the charge on the superconducting island in the middle of 
the hexagon. 

In the experiment the differential resistance in the fiux-fiow regime has been measured as 
a function of the gate voltage. When fixing the current through the array, clear oscillations 
in the differential resistance are observed. The measured period, however, corresponds to 
e, half the value expected from theory. The factor of two arises from tunneling of quasi- 
particles which effectively limit the quantum phase difference to values in the range — 7r/2 
to 7t/2. (Quasi-particle tunneling changes the vortex phase difference by tt and becomes 
favorable as soon as the induced charge equals e/2, i.e. when the phase difference equals 
vr/2.) 

It is possible to derive the AC effect from the Quantum Phase Model. In a more transpar- 
ent way, although less rigorous, one starts from the representation of the partition function 
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as a path integral over the phases and charges (see Eg. (P2|) ) with the inclusion of a uniform 
background charge, i.e. q — > q — q^- The term in the action of the QPM which is relevant 
for the AC effect is the one which is linear in 0j in Eq. ( |D2| ). If a vortex is present, the 
configuration of the phases will be related to its position r(r) and, by going over the 



same steps that lead to the single vortex action in Eq. |7T|, there is an additional term to the 
effective action equal to 



Sac = -i 



E^^'^^r.Ve[r,-r(r)] (80) 



Eq.(pOD defines a pseudo-charge gauge field for the vortex 

2eAQ(r) = q,^z x r/r^ 

which is singular at the origin of the vortex. Thus the external charge act like a vector 
potential for the vortex. The phase factor x implied by Eq.(|80D is 

X = 2vrE|^ (81) 

icr 

where the sum extends to all islands enclosed by the trajectory T . 



3. Block oscillations 



Electrons in metals move in the periodic potential created by the positively charged ions. 
The electron wave functions overlap and energy bands are formed. A constant electric field 
accelerates electrons, but in the absence of scattering, electrons would be Bragg reflected 
at the Brillouin zone edges. Electrons then undergo an oscillatory motion in space (Bloch 
oscillations). No charge would be transported. In metals scattering takes place before 
the electrons can reach the zone edge so that Bloch oscillations do not appear. In semi- 
conductor superlattices ||279|| Bloch oscillations have been observed because of the larger 
superlattice period and because of less scattering in the controlled fabricated structures. 
Coherent Cooper pair tunneling in current bias Josephson junctions leads to a phenomenon 
analogous to Bloch oscillations [p80| . 

Vortices in a periodic potential should also from energy bands. It is possible to study 
Bloch oscillation for vortices ||267|| in a quasi- ID Josephson array that is a few cells wide 
and 1000 cells long. A sketch of the sample layout is shown in Fig.|36|. For low densities, 
the bus-bars force the vortices to move in the middle row so that they experience a purely 
ID sine potential. For a free vortex the energy depends quadratic on the wave vector 
k: E{k) = k'^/2My, which equals E{k) = 2Ec at the Brillouin zone edge. In a periodic 
potential, energy gaps open up at the zone edges. The gap is equal to the Fourier coefficient 

For a sine potential ^•jEj sm{2'7ix), the gap is then •yEj. 



of the lattice potential 
Thus, vortices in arrays form energy bands with a bandwidth of the order Eq and an 



energy gap of ^■yEj as illustrated in Fig.^. Assuming the lowest band to be cosine-shaped 
{E{k) = hEc{l — cos(A;))), the equation of motion F = hdk/dt is: 



^dk 
dt 



$0-^ ~ Tiv{k) where 



^sin(/.). 



(82) 
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As defined before, / denotes the applied current per junction, t] the phenomenological vis- 
cosity and v{k) the average vortex velocity. 

In the absence of damping, with a small current applied, the wave vector changes linearly 
in time: the vortex thus reaches the Brillouin zone edge where it will be Bragg reflected. 
This Bragg reflection results in an oscillatory motion in fc-space. On average the vortex 
velocity is zero and the time it takes the vortex to complete one oscillation follows from 
At = Ak/ < dk/dt > with Ak = 2tt. The corresponding Bloch oscillation frequency (ub) 
is: 

= ^ , (83) 

and the amplitude of the oscillation is 

^ = 1^ = ^h- (84) 

J (Pol Ej Til ^ ' 

When biasing an array with 1000 junctions with currents of the order of yuA, Bloch fre- 
quencies are in the range 1-10 GHz. Since Ec ~ Ej and Id I is typically 100, the Bloch 
oscillations extend over 10 cells. 

A characteristic feature of Bloch oscillating vortices is a nose-shaped form of the dc 
current- voltage characteristic (see Fig.^8|). For very small bias, there is a small supercurrent 
because vortices need to overcome the energy barriers near the array edges (finite-size effect). 
Just above the depinning current, any amount of dissipation prevents vortices from reaching 
the zone edges. Bloch oscillations do not exist in this regime and an increase of the current 
yields an increase of the measured voltage across the array. 

When increasing the current beyond some point, dissipation is not strong enough to 
prevent the vortices from reaching the zone edges. Bloch oscillations are now possible. In 
the I-V characteristic a sudden decrease of the voltage is then expected with a negative 
differential resistance: the oscillating vortices do not contribute to the net transport of 
vortices through the array. Eqs.(p2D can be solved with the result 

V{I) =7171— ^ if / < Jo (85) 

e /q 



V(I)=mi— 

e ^ 



where Jq = rjEcW/ (2$o) and n is the one-dimensional vortex density (=number of vortices 
divided by the number of cells in the direction of motion). The solid line in Fig.|3^ is a fit to 
these equations. (The line is offset by a small positive current to correct for the depinning 
current). Although the overall shape of the experimental curve resembles that of the theo- 
retical prediction, the experimental value of the band width {Ec is the previous discussion) 
is one order of magnitude smaller than expected. This discrepancy is not understood, but 
is consistent with a vortex mass that is larger than the calculated, quasi-static vortex mass. 
Note that the data extracted from the study of the macroscopic quantum tunneling also 
indicated the same trend for the experimental vortex mass. 



56 



Additional information on Bloch vortices can be obtained by irradiating the sample with 
a microwave signal. Steps occur in the I-V characteristics when the external frequency 
locks to the Bloch frequency. Surprisingly, the experiments show that the Bloch frequency 
depends on the vortex density. This dependence is not accounted for by the independent- 
vortex model presented above (see Eq. (p3D) and suggests a collective oscillation of the 
ID vortex chain (see next section). We will come back to this issue when discussing the 
formation of a Mott insulator in these quasi-lD arrays. 



E. One-dimensional vortex localization 



In Section [II G , we already introduced quasi-one dimensional arrays as model systems 



for the study of interacting bosons in one dimensions. For an ideal periodic potential and in 
the absence of interactions between the particles, the solution of the Schrodinger equation 
consists of Bloch waves that extend throughout the whole chain. 

In quasi-lD Josephson arrays, however, vortex localization can occur in several ways. 
We already discussed the existence of Bloch oscillations when an external force is exerted on 
a vortex. One can view this oscillatory motion in space as localization of individual bosons 
analogous to Wannier-Stark localization of electrons ||282||: the extend of the wave function 



is decreased when the external force on the quantum particle is increased. 

In the next two subsections, we treat two other cases where boson localization occurs. 



Commensurability effects ||283|| may lead to localization of vortices in quasi- ID arrays. The 
repelling interaction between vortices plays a crucial role in this so-called Mott-localization 
PPI . Another mechanism to localize Bloch waves is disorder and this phenomenon is known 



as Anderson localization [ P84|| . Two important remarks should be made at this ponit. First, 
vortex localization, as discussed in the next subsection, has been studied by measuring 
the zero-bias resistance. In this case only a very small current is applied in contrast to 
the experiment showing the Bloch oscillations. Second, one should realize that in contrast 
to localization in electronic systems, vortex localization leads to a zero resistive state. In 
superconductors motion of vortices is the cause of dissipation. If they are localized, the 
sample is superconducting. 



1. Mott insulator of vortices 

More than ten years ago, localization of bosonic particles with a short-range repulsive 
interaction has been studied theoretically by Fisher et al. [^. They found a Mott insulating 
phase for commensurate filling and a superfluid phase for incommensurate filling. Strong 
disorder destroys the Mott phase and there is the possibility to have a Bose glass (for 
theoretical studies of the phase-diagram of bosons on a chain see Refs. ||118| , |119| , |285|] and 
references therein). 

Experimentally, Mott localization of vortices has been studied by van Oudenaarden et al. 
||263| , |265|| . They explored the influence of the interaction strength, the bandwidth, the sample 
geometry and temperature on the stability of the Mott states. First, we will summarize the 
experimental results. Then, we discuss their experiment in the context of a recent theory of 
Bruder et al. [|8|]. 
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The experiment consists in measuring the zero-bias resistance vs. magnetic field for 
quasi-one dimensional arrays of different lengths and Ej/Ec ratios. It is convenient to 
define a ID frustration n, which is associated with magnetic field piercing through a cell of 
area W (in units of the lattice constant): n = WB/^q = Wf. In Fig.^ the results are 
shown for arrays of three and seven cells wide, i.e., W = 3 and 7 respectively. The plot 
for the W = 7-sample is mirrored with respect to the x-axis for clarity. For both samples, 
clear dips occur for certain values of the 2D frustration index /, i.e., for certain values of 
the vortex density. When plotted vs. this ID frustration index, the ID nature of the sharp 
dips becomes visible. Dips are found at the same rational values of n = 1/3, 1/2, 1, 2. 

More detailed measurements of the resistance dips show that they are not infinitely 
sharp. There is a certain window of n in which the resistance vanishes and the vortex 
chain is pinned. The interaction energy proportional to Ej (see below) dominates the 
bandwidth proportional to Eq in this regime. Beyond this window the resistance increases 
sharply, indicating that the vortex chain is depinned and that the bandwidth dominates the 
interaction energy. From this consideration, one expects the window to be larger for samples 
with a larger Ej/Ec ratios. In the experiment, this dependence has indeed been observed. 

Around commensurate filling, the system is incompressible: small changes of the mag- 
netic field (i.e. n) do not lead to a change in the number of vortices in the chain. This process 
costs a certain energy, called the Mott gap. By analyzing thermally activated transport in 
the Mott states, the value of this gap can be deduced from the experiments and values 
in the order of Kelvins are reported. The infiuence of the array length was also studied. 
No significant differences were observed in arrays with lengths larger than 200 cells. This 
observation demonstrates that edge effects do not play an important role and that the long 
arrays are indeed one- dimensional systems. 

A quantitative theory of the commensurate-incommensurate transition of a vortex chain 
in quasi-one-dimensional Josephson arrays has been formulated by Bruder at al. |p86|| . They 
showed that the transition to the incommensurate state is due to the proliferation of soliton 
excitations of the vortex chain. Since the range of the interaction between the vortices is 
much longer than the inter-vortex distance, solitons consist of many vortices, and possess a 
large effective mass. The transfer of one flux quantum between the array edges is then due 
to soliton propagation through the sample. The number of solitons necessary to transfer 
one vortex is equal to the ratio of the periods of the vortex lattice and the junction array. 

The analysis of Bruder at al. ||286|| focuses on determining the energy barrier Er of 
the observed thermal activation in the resistance vs. temperature curves. Although this 
approach is of semi-classical nature, quantum effects are crucial in relating the parameters 
of the effective theory (expressed in terms of soliton excitations) to the microscopic couplings 
{Ec and Ej) of the Josephson array. 

In the commensurate phase, there are two contributions to Er. The flrst contribution 
comes from the activation energy of a soliton and the second term summarizes the boundary 
pinning energies. This boundary effect can be understood as follows: Because of commen- 
surability, the process of vortex flow through the array can be viewed as the motion of a 
rigid vortex chain. Therefore, the vortex chain cannot adjust itself to the boundary pinning 
potential. The potentials produced by the two array ends both contribute to Er and the 
relative phase of these two contributions depends on whether the total flux piercing the 
junction array equals an integer number of flux quanta or not. Consequently, the second 
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-boundary pinning- term to Er oscillates with the magnetic flux piercing the array. In Fig.|0 
Er is plotted for the case of small boundary pinning while the opposite situation is shown 
in Fig.p|. The short period oscillations are determined by boundary-pinning term while the 
vanishing of Er at the edge of Mott lobe is driven by the soliton energy. 

In the incommensurate state, the vortex chain is compressible, and can adjust itself to 
the boundaries of the array. As a result, the main contribution to the activation energy is 
due to the boundary pinning potential and the elastic energy. 

The comparison between theory and experiment is shown in Fig.^. The theoretical 
results lead to an estimate for the soliton energy much larger than the boundary pinning 
which is of the order of ~ 0.57^. The theoretical value of the activation energy 

Er ~ 8K 

is in very good agreement with the experiments. 

Two observations provide additional support for the theory of Bruder et al. First, the 
regions of n corresponding to the Mott phase are extremely narrow suggesting, at least in 
the conventional non-interacting picture, a weak interaction between the particles. Conse- 
quently, within the Mott phase the activation energies for particle transport are expected to 
be small as well. However, the observed value of Er is one order of magnitude larger than 
the energies Ec and Ej, which determine the single- vortex band spectrum. A second feature 
in favor is the strong oscillating behavior of the resistance (with a period proportional to 
the inverse length of the chain) outside the Mott region. These oscillations would not be 
expected in a model of almost-free quasi-particles within the delocalized phase. 



2. Anderson localization of vortices 

In 1958, Anderson [ ^84j| showed that disorder has a dramatic influence on transport 
properties. Disorder reduces the spatial extent of wave functions to such an extent that 



transport can completely be blocked. Three years later Mott and Twose ||287|] showed that 
ID systems are particularly susceptible for disorder: even weak disorder leads to strong 
localization. Many studies on Anderson localization have been performed on 3D and 2D 
samples. One-dimensional model systems are harder to find. Josephson-junction arrays have 
the great advantage that disorder can be introduced in a controlled way. Experimentally, 
Anderson localization of vortices has been studied by the Delft group |264] | . Their results 
will be outlined in this subsection. 

In the experiment by the Delft group disorder has been introduced by constructing 
superlattice structures. The superlattice is formed by replacing all the junctions of a column 
by junctions that are twice as large. Consequently, these barrier junctions have a Josephson 
energy that is two times larger than that of the adjacent junctions. The barrier junctions 
yield a peak in the potential landscape for vortices traveling through the array. Numerical 
calculations show that this barrier is 1.7Ej, which is about one order of magnitude larger 
than the energy barrier for cell-to-cell motion. 

A perfect superlattice structure is made by introducing columns of barrier junctions on 
a distance of exactly 10 lattice cells. For an array of length 1000, this means that there are 
100 columns that have been changed. Disorder is now introduced by changing the distance 



59 



between two barriers Samples with different amounts of disorder have been fabricated. In 
the least disordered samples (labeled with 6 = 1) the barriers were separated by 9, 10 or 
11 lattice cells with equal probabilities. In other disordered samples (labeled with 6 = 2) 
barriers placed at distances 8, 9, 10, 11, or 12 lattice constants again with equal probability. 
All samples contained 100 barriers. 

The vortex quantum properties are probed by measuring the zero-bias resistance as a 
function of temperature for the perfect periodic array as well as for the disordered arrays. 
Since the topic of interest is the study of quantum transport, the vortex density needs to have 
a non-commensurate value to avoid the Mott state as discussed in the previous subsection. 
The result for n = 0.44 is shown in Fig.^. At high temperatures all three arrays show the 
same behavior: transport is thermally activated with an energy barrier of 3Ej, a factor of 
two larger than the expected value. When the temperature is lowered, however, a significant 
difference is observed between the periodic sample and the two disordered samples. For the 
perfect periodic sample a finite resistance is measured at the lowest temperatures. In this 
regime the resistance is independent on temperature indicating vortex transport by quantum 
tunneling. In contrast, the resistances of the disordered samples have dropped below the 
measuring accuracy of the set-up. Thus, in the perfect array vortices are mobile whereas 
they are localized in the disordered arrays. 

The zero-bias resistance has been studied for several values of the vortex density. For 
n < 0.3, the zero-bias resistance is too small to be resolved for the periodic array. In the 
range 0.3 < n < 0.8, the resistance of the periodic array is significantly larger than that 
of the disordered arrays and this is the region where vortex localization occurs. But for 
even larger vortex densities the behavior of all three arrays is almost the same showing 
a flattening of the resistance at the lowest temperatures. In all three arrays vortices are 
now mobile. At these high vortex densities the distance between them is small and their 
repulsive interaction can no longer be neglected. The experiment shows that in this case 
delocalization occurs. 

Above we have discussed the experiment in terms of Anderson localization which strictly 
speaking occurs when the interaction between the bosons is very weak; i.e., when the bosons 
act as independent particles that are localized for arbitrarily weak disorder. In the exper- 
iment, vortex densities are large. Disorder now competes with the interaction strength. 
A sufficiently strong interaction can delocalize the particles, whereas strong disorder will 
localize them again in a Bose-glass phase. To distinguish between Anderson localization 
and localization in a Rose glass, more measurements are needed. Experiments should be 
performed on arrays with fixed disorder but different Ej to clarify the role of the interaction 
strength. 

IV. FUTURE DIRECTIONS 

In this last chapter of the review, we discuss some future directions for research on 
Josephson networks. Of course, additional information on the quantum nature of Joseph- 
son networks can be obtained by using new measuring techniques (e.g. the ac-measuring 
method that has successfully been applied to classical arrays) and better samples (e.g. ar- 
rays with a well-controlled environment). Here, we outline the concepts behind two experi- 
ments -persistent vortex currents and vortex quantum Hall effect- of which some theoretical 
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calculations are available. We summarize the main ideas as well as the experimental re- 
quirements/improvements for observation of these effects. Future lines of research may also 
involve the use of Josephson networks as model systems in unexplored areas of physics. Bio- 
physics may be such a field and an interesting example in this respect is vortex transport in 
ratchet arrays. Undoubtedly, the most exciting new line of future research is quantum com- 
putation. In this paper, we do not have space to treat quantum computation with Josephson 
circuits in great depth. We therefore only present the concept and summarize the current 
status of art. 



A. Persistent vortex currents 

Persistent currents in rings made of normal metals are a manifestation of the Aharonov- 
Bohm effect for quantum coherent electrons in a multiply connected geometry. In the absence 
of any driving current, the flux through the ring induces the motion of charge carriers that 



can be detected at low temperatures ||288|| . Charge-fiux duality indicates that in Josephson 
rings a persistent current of vortices may be expected generated by the charge induced on 
the inner island. A persistent vortex current leads to a persistent voltage across the ring 
and as such would be manifestation of the Aharonov-Casher effect discussed in Sec. [ill D| . 
Although vortex interference in an open Josephson circuit has been reported by Elion et 



al. there are no experiments reporting the existence of persistent vortex currents in 



Josephson Corbino circuits. The proposed setup to measure persistent vortex current is 



shown in Fig. 44 



When the mean free vortex path is long enough, we can neglect the dissipation in the 
dynamics and the Hamiltonian of a vortex in a discrete Josephson ring (with N junctions) 
can be written as 

1 fr. *0. 



H = P -Q^ (87) 



where P is the vortex momentum and the charge on the inner island of the ring. The 
vortex dynamics is quantized and the set of discrete energy values is given by 



where i> is an integer and where Ce// agrees with the continuum result of Ref. ||289|] and with 



Ckin in Ref. ||278|| . Thus, a Josephson ring with a vortex trapped inside acts like a perfect 
capacitor. A similar situation arises in a single classical junction in the absence of screening. 
With exactly one vortex trapped, its critical current is zero. The difference lies in the value 
of the effective capacitance, which in quantum rings differs from the geometrical one. For 
example, for a ring consisting of a square 2D array Ce// = {N'^/2)C and for a ID ring with 
Aj < N, Ceff = {2N/7r^Aj)NC. 

Although this change in capacitance can in principle be measured, quantum rings should 
exhibit a more interesting phenomenon: persistent vortex currents. This can easily been seen 
from the Hamiltonian of Eq. (P7|) since it has the same from as the Hamiltonian for electrons 



in a metal ring. Duality arguments then indicate the existence of a persistent voltage due 
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to the persistent motion ||225|,|277|,|290|| of a vortex induced by Qx- The basic reasoning is 



as follows: With some disorder, small gaps open up in the energy spectrum of Eq.(88) and 
energy bands form. With no current applied, the charge dictates the quantum dynamics 
of the ring. It determines the vortex velocity and therefore the voltage across the ring 
because the vortex speed is proportional to dEn/dQx- As a result, a persistent voltage 
appears across the ring which is periodic in Qx with period 2e. 

Within the free-particle model, the maximum voltage can be estimated. The voltage 
V due to a circling vortex with velocity v is equal to ^qv/N. The velocity is equal to 
{N/(^Q)dEy/dQx- Then for = and = e, the maximum voltage equals e/Ceff = 
e/M„ (<l>o/A^)^. For a ring consisting of a square 2D array V^ax = {2e/CN'^). With = 10, 

C = 1 fF, Vmax ~ 3 /iV. 

In the picture presented above, the vortex is treated as a free, non-interacting particle. 
A detailed description of vortex dynamics, however, is only possible if one considers its 
dissipative environment. As discussed before, an important dephasing mechanism is the 
existence of linear spin-waves. This coupling leads to damping and hence to a finite phase- 



coherence length for vortices. This problem was studied in Ref. |p91|| . For the persistent 
voltage one obtains 



.2 V) = E^=i^sin(27rz/g./2e)exp(-^,) 

^ ' /5 l + 2E^=iCos(27rz/Qj2e)exp(-^,) ^ ^ 



where 5*,^ corresponds to the action in the sector of u winding numbers. In the adiabatic 



mass approximation in the limits Ej ^ Eq and -Rjv ~ Rk, the action reduces to that of a 
free particle with mass M„ = n'^/AEc- The saddle point action is 

S'. = ^{^Nr (90) 

and the persistent voltage is a sawtooth function of at zero temperature and a sine- 
like function at higher temperatures. It depends on the system size through the ratio of the 
radius N/ n to the thermal wavelength of the vortex At = \J2tiJm^ = ^SPEc/tc- 

To go beyond the adiabatic mass approximation, the actions can be calculated with 
the full non-local kernel (see Eq. |7l|) that incorporates the effect of inelastic processes due 
to the interaction with spin-waves. An important consequence of the quantum corrections 
is that the vortex dephasing length 

= oOp^/T (91) 

increasing with the ratio Ec ~ Ej (^ gets larger on approaching the S-I transition). For 
Ec Ej IK, L^p may be as large as 100 lattice constants at T= 10 mK, which is much 
larger than the thermal wavelength Ay for the vortex in the adiabatic mass approximation. 
In the limit ^ <^ and <^ oOp the adiabatic result is recovered. If ^ ~ 1 a larger 
persistent voltage is found. If the ratio Ej/Ec is reduced, the coherence length ^ grows 
beyond the radius of the system, and the persistent voltage grows even more. As compared 
to the adiabatic mass limit, the persistent voltage including the vortex-spin- wave coupling 
is always larger. For Ej= Ec=l K, T=10 mK and = 10, a persistent voltage in the 
microvolt range is expected, which is observable. 



62 



Other interesting issues involve the effect of the underlying lattice and of disorder on the 
persistent voltage. They both induce backscattering and the opening of gaps in the band 
structure E{Qx). One expects Zener tunneling across the band gaps to occur if the voltage 
Vx is switched on fast, yielding a higher transient current that relaxes to the persistent 
current after some relaxation time. 

Experimentally, fabrication of quantum rings is difficult, but does not seem to be impos- 
sible. First of all vortices should have long mean free paths. As we have seen, 2D arrays 
exhibit only a small window for which free propagation of vortices is possible. Purely ID 
discrete Josephson rings with one vortex trapped | |228| j seem to be more promising in this 
respect since for Aj > 1 there is no energy barrier for vortex motion. Coupling to spin-waves 
can be small and long mean free vortex paths are expected [p92|| . Secondly, Josephson rings 
should be decoupled from their environment. This can be achieved by placing high-Ohmic 
resistors or alternatively arrays of junctions in the leads close by. 

A third restriction comes from the charging energy. Calculations on the continuous 
Josephson system ||293|| show that the temperature has to be smaller than e^/(27r^Ce/j) in 
order to observe quantum vortex dynamics. Temperatures of the order of 100 mK therefore 
require Ce// to be 1 fF. This requirement indicates that aluminum junctions have to be 
smaller than 0.01 yum^ and that the whole ring structure can not be made too large. The 
capacitances to ground would otherwise be too dominant. The fourth restriction comes from 
the shadow-evaporation method itself. Both the wire connecting the central island in the 
middle of the ring and the gate capacitor -preferable situated underneath or on top of the 
central island- have to be made in separate fabrication steps. Alignment and good electrical 
isolation between the different layers have to be established. 



B. The Quantum Hall effect 



Quantum electron transport in two-dimensional systems in the presence of an applied 
magnetic field is one of the most intensively investigated areas in condensed matter. When 
the filling factor (ratio between the electron density and the density of flux quanta ) is of the 
order of one, the Quantum Hall effect (integer and fractional) occurs [ p94|| . In the previous 
sections we showed that charges and vortices behave as quantum particles hopping coherently 
in the artificial two-dimensional space created by the Josephson array. The duality between 
charges and vortices can be shown also in the presence of external frustration; magnetic field 

3|- |297| | it was suggested 



for charges and offset charges for vortices. In a series of papers 
that a Quantum Hall effect could be observed in Josephson arrays both for charges and 
vortices. The three proposals address different regimes: Nazarov and Odintsov p95| describe 



the possibility of Hall states for Cooper pairs while the authors in Refs. [P96| , P97]| consider 
the Hall fluid of vortices. 

In the limit Eq » Ej and in the case of very low density, charges behaves as a dilute 
Bose gas with strong repulsion. Under these conditions, analytic and numerical calcula- 
tions |p95|] support the idea that in a magnetic field Cooper pairs form Laughlin-type in- 
compressible states (a Cooper pair fluid): the charge density changes in a stepwise function 
by changing the external parameters. The incompressible states give rise to the quantization 
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of the Hall conductance 



'xy 



h 



where the filling factor u = q/f is given by the ratio between the charge density q and the 
magnetic frustration /. According to Odintsov and Nazarov two sets of Hall plateaus exist. 
One corresponds to the fractional quantum Hall effect with u = 2m (m integer). The other 
corresponds to the integer quantum Hall effect with v = 1/2 {I integer). 

The opposite limit in which the Josephson energy dominates, has been considered in 
Refs. [|296| , p97|| . In this case vortices condense to form a quantum Hall fluid and the transverse 
conductivity is now given by: 

a^y = 2m— . 

Despite the similarities highlighted here, there are important differences as compared to 
the "electronic" case. Both charges and vortices are bosons, moreover they are interacting 
particles. Interactions modify the results. For instance, in the vortex case the logarithmic 
interaction changes the longitudinal response [ 297)1 as well. 

Up to now there is no experimental evidence for the quantum Hall effect in Josephson 
arrays. For the Cooper-pair fluid, the random offset charges form a serious obstacle for 
observation of the Quantum Hall states. For the vortex fluid the situation is less clear. 
Additional theoretical work is required to locate the region in parameter space where the 
Hall fluid is the ground state (as opposed to the Abrikosov lattice for example). Effects 
related to disorder, quantum correction of the mass, and dissipation should all be taken into 
account. 



C. Quantum Computation with Josephson junctions 

Quantum Computation (QC) has recently excited many scientists from various different 
areas of physics, mathematics and computer science. In contrast to its classical counterpart, 
quantum information processing is based on the controlled unitary evolution of quantum 
mechanical systems. The great interest in this fleld is certainly related to the fact that 
some problems which are intractable with classical algorithms can be solved much faster 
with QC. Factorization of large numbers as proposed by Shor is probably the best known 
example in this respect. This section briefly reviews the recent work in this fleld using 



Josephson nano-circuits. For excellent reviews devoted to QC we refer to Refs. ||298| , |299 
Furthermore, many elementary books on quantum mechanics treat the physics of two-level 
systems. A review devoted to the implementation of quantum computation by means of 
Josephson nano-circuits just appeared [pOO|| . 



The elementary unit of any quantum information process is the qubit. The two values 
of the classical bit are replaced by the ground state (|0 >) and excited (|1 >) state of a 
two-level system. (Note that it common to adopt the spin-1/2 language as we will do here.) 
Already at this stage a fundamental difference between classical and quantum bits emerges. 
While information is stored either in or in 1 in a classical bit, any state \ip{t) >= a{t)\0 > 
+bit)\l > can be used as a qubit. 

Manipulations of spin systems have been widely studied and nowadays NMR physicists 
can prepare the spin system in any state and let it evolve to any other state. Controlled 
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evolution between the two degenerate states |0 > and |1 > is obtained by applying resonant 
microwaves to the system but state control can also be achieved with a fast DC pulse of high 
amplitude. By choosing the appropriate pulse widths, the NOT operation can be established 

|0 > — ^ |1 > (92) 

11 > — ^ 10 > 



or the Hadamard transformation 

|0 > ^ -L(|0 > +|1 >) (93) 

11 > ^^(|0>-|1». 

These unitary operations alone do not make a quantum computer yet. Together with one- 
bit operations it is of fundamental importance to perform two-bit quantum operations; i.e., 
to control the unitary evolution of entangled states. Thus, a universal quantum computer 
needs both one and two-qubit gates (it has been shown that most of two-qubit gates are 
universal [p01|| ). One example of a two-qubit gate is the Control-NOT operation: 



00 > - 


|00 > 


01 > - 


|01 > 


10 > - 


111 > 


11 > - 


|10 > 



(94) 



The unitary single-bit operations and this Control-NOT operation are sufficient for per- 
forming all tasks of a quantum computer. Therefore, quantum computers can be viewed 
as programmable quantum interferometers. Initially prepared in a superposition of all the 
possible input states, the computation evolves in parallel along all its possible paths, which 
interfere constructively towards the desired output state. It is this intrinsic parallelism in the 
evolution of quantum systems that allows for exponentially more efficient ways of performing 
computation. 

It is of crucial importance that qubits are protected from the environment, i.e., from 
any source that could cause decoherence [P02|| . This is a very difficult task because at the 
same time one also has to control the evolution of the qubits, which inevitably means that 
the qubit is coupled to the environment. In quantum optics experiments, single atoms are 
manipulated which are almost decoupled from the outside world. Large-scale integration 
(needed to make a quantum computer useful) seems to be on the other hand impossible. 
Qubits made out of solid-state devices (spins in quantum dots or superconducting nano- 
devices), may offer a great advantage in this respect because fabrication techniques allow 
for scalability to a large number of coupled qubits. 

At present different proposals have been put forward to use superconducting nano- 
circuits ||303H307| for the implementation of quantum algorithms. Depending on the op- 
erating regime, they are commonly referred to as charge [|303| , |304| , ^07| , ^08|| and flux [ ^05|j306|| 
qubits. We briefly discuss both approaches and summarize the experimental advances made 
so far. 
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Charge qubits |P03|J304|| - In this case the qubit is reahzed by the two nearly degenerate 



charge states of a single electron box as shown in Fig.|45|. They represent the states |0), |1) of 
the qubit. In the computational Hilbert space the ideal evolution of the system is governed 
by the Hamiltonian 

H = -AEeh,i((|0 >< 0| - |1 >< 1|) - ^(|0 X 1| + |1 X 0|) (95) 

where AEch,i = E^hiux — 1/2). Any one bit operation can be realized by varying the 
external charge n^. and, in the proposal of Ref. |^03[ by varying the Josephson coupling as 



well. Modulation of Ej is achieved by placing the Cooper-pair box in SQUID geometry. 
The advantage of this choice is that during idle times the Hamiltonian can be "switched 
off' completely eliminating any trivial phase accumulation which should be subtracted for 
computational purposes. 

As discussed before, a quantum computer can be realized once two bit gates are imple- 
mented. The Karlsruhe group has proposed an inductive coupling between qubits which 
lead to a coupling of the type 

i/c = -i^L^'VJ^) . (96) 

This type of coupling is very close in spirit to the coupling used in the ion-trap implemen- 
tation of QC. The main advantage of this choice is that qubits are coupled via an infinite 
range coupling and that the two bits can easily be isolated. A different scheme has been 
proposed in Ref. ||304|| . They emphasize the adiabatic aspect of conditional dynamics and 



suggest to use capacitive coupling between gates as to reduce unwanted transitions to higher 
charge states. The coupling reads He = —Eccr^^^cr^p and the qubit is now defined as a finite 
one-dimensional array of junctions. By means of gate voltages applied at different places in 
the array the bit-bit coupling can be modulated in time and a control-NOT can be realized. 
The experiments on the superposition of charge states in Josephson junctions [p09|j310 



and the recent achievements in controlling the coherent evolution of quantum states in a 



Cooper pair box ||311|| render superconducting nanocircuits interesting candidates to imple- 



ment solid state quantum computers. The experiment by Nakamura et al. goes as 

follows. Initially, the system is prepared in the ground state. Appropriate voltage pulses 
bring the system in resonance so that the two charge states are in a coherent superposition 
a(i)|0 > +h{t)\l >. The final state is measured by detecting a tunneling current through 
an additional probe-junction. For example, zero tunneling current implies that the system 
ended up in the |0 > state, whereas a maximum current indicates that the final state cor- 
responds to the excited one. In the experiment the tunneling current shows an oscillating 
behavior as a function of pulse length, thereby demonstrating the evolution of a coherent 
quantum state in the time domain. 

Nakamura et al. also estimate the dephasing time and report it to be of the order of 
few nanoseconds. The probe junction and 1/f noise presumably due the motion to trapped 
charges are the main source of decoherence. In their absence, the main dephasing mech- 
anism is thought to be spontaneous photon emission to the electromagnetic environment. 
Decoherence times of the order of 1 /is should then be possible. 



Phase qubits ||305|J306|| - A qubit can also be realized with superconducting nano-circuits 



in the opposite limit Ej ^ Eq- An rf-SQUID (a superconducting loop interrupted by a 
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Josephson junction) provides the prototype of such a device. The Hamiltonian of this system 
reads 



Here, L is the self- inductance of the loop and the phase difference across the junction 
(27r$/$o) is related the flux $ in the loop. The externally applied flux is denoted by The 
charge Q is canonically conjugated to the flux $. In the limit in which the self-inductance 
is large, the two first terms in the Hamiltonian form a double-well potential near $ = $o/2. 
Also in this case the Hamiltonian can be reduced to that of a two-state system. The term 
proportional to measure the asymmetry of the double well potential and the off-diagonal 
matrix elements depend on the tunneling amplitude between the wells. By controlling the 
applied magnetic field, all elementary unitary operations can be performed. 

In order to fulfill various operational requirements more refined designs should be used. 
In the proposal of Mooij et al. [p05|] , qubits are formed by three junctions (see Fig. ^6]). Flux 
qubits are coupled by means of flux transformers which provide inductive coupling between 
them. Any loop of one qubit can be coupled to any loop of the other, but to turn off this 
coupling, one would need to have an ideal switch in the flux transformer. This switch is to be 
controlled by high-frequency pulses and the related external circuit can lead to decoherence 
effects. 



At present, both the Stony Brook group (Friedman et al. ||312|| ) and the Delft group 
(Van der Wal et al. [ pl3|] ) have demonstrated superposition of two magnetic flux states 
in superconducting loops. One state corresponds to the magnetic moment of /lA-currents 
flowing clockwise whereas the other corresponds the same moment but of opposite sign 
due to the current flowing anti-clockwise. Coherent quantum oscillations have not yet been 
detected. To probe the time evolution, pulsed microwaves instead of continuous ones have to 
be applied. Observation of such oscillations would imply the demonstration of macroscopic 
quantum coherence (MQC). It is called macroscopic because the currents are built of billions 
of electrons coherently circulating within the superconducting ring. 

There are two main differences between the approaches of the two groups. The Stony 
Brook group uses the excited states of an RF SQUID. The Delft circuit consists of a three- 
junction system and the continuous microwaves induce transitions between the ground state 
and the first excited state only. The three-j unction geometry has the advantage that it can 
be made much smaller so that it is less sensitive to noise introduced by inductive coupling 
to the environment. Nevertheless, recent insights indicate that in all designs put forward 
so far the measuring equipment destroys quantum coherence. The meter is believed to be 
main obstacle to study the 'intrinsic' decoherence times. Future work must evaluate the role 
of the measuring equipment and new measuring schemes should be developed in order to 
study MQC in Josephson loops. 
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APPENDIX A: ARRAY FABRICATION AND EXPERIMENTAL DETAILS 

There are two types of junctions arrays: proximity coupled arrays and arrays made of 
Josephson tunne/ junctions. The proximity coupled arrays consist of superconducting islands 
(Nb or Pb) on top of a normal metal film (Cu) and are solely used for the study of classical 
phenomena. The main reason for this is the low (less than 1 Ohm) normal-state resistance 
of the junctions. 

At present there are two technologies to fabricate arrays with Josephson tunnel junc- 
tions. Commercial niobium junctions are fabricated using a trilayer with aluminum oxide as 
insulating barrier. Reliable niobium junctions, however, are too large to observe the quan- 
tum effects discussed in this review. Quantum arrays are built up of all aluminum, tunnel 
junctions. These Josephson junctions are fabricated with a shadow-evaporation technique 
| 314|| . We will now outline the most recent fabrication technique as used in the Delft group. 

Samples are fabricated on silicon substrates with an insulating Si02 top layer. In the first 
step two resist layers are spun on the substrate. The lower layer is a solution of PMMA/MAA 
copolymer in acetic acid and the upper layer is a solution of PMMA in chlorobenzene. The 
resist sandwich is baked at 180 °C for one hour. Then the sample mask is written by 
high-resolution electron-beam lithography at 100 kV. After writing, the exposed resist is 
developed in a 1:3 mixture of MIBK and 2-propanol for one minute. The solubility of the 
lower resist is larger than the upper layer, which leads to an undercut in the lower resist 
layer. This undercut is necessary for the formation of free-hanging bridges below which the 
junctions are formed in the shadow evaporation technique. 

After mask definition, a 24 nm thick aluminum layer is evaporated under a given angle. 
Then the aluminum layer is exposed to pure oxygen at a controlled pressure. By changing 
the pressure the thickness of the aluminum oxide barrier is varied. In the second evaporation 
step (the sample is not taken out of the vacuum) a 40 nm thick aluminum layer is evaporated 
under the opposite angle. After this step, the tunnel junction is formed. The remaining resist 
layers with the unwanted aluminum on top are removed by rinsing the sample in acetone. 

We end this appendix with some remarks on the measuring set-up. Arrays are measured 
in a dilution refrigerator inside /i-metal and lead magnetic shields at temperatures down to 
10 mK. To protect the arrays from high energy photons generated by room-temperature noise 
and radiation, extensive filtering and placing the arrays in a closed copper box are minimum 
requirements. Therefore, a typical set-up for the measurements of quantum arrays has the 
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following characteristics. At the entrance of the cryostat, electrical leads are filtered with 
radio-frequency interference (RFI) feedthrough filters. Arrays are placed inside a closed, 
grounded copper box (microwave-tight). All leads leaving this box are filtered with RC 
filters for low-frequency filtering {R = 1 kfl and C = 470 pF) and with microwave filters. A 
microwave filter consists of a coiled manganin wire (length ~ 5 m), put inside a grounded 
copper tube that is filled with copper powder (grains < 30 /im). The resistance of the wire 
in combination with the capacitance to ground via the copper grains provide an attenuation 
over 150 dB at frequencies higher than 1 MHz. The copper box with the RC and microwave 
filters is situated in the inner vacuum chamber and is mounted on the mixing chamber in 
good thermal contact. 



APPENDIX B: TRIANGULAR ARRAYS AND GEOMETRICAL FACTORS 

In comparing properties of square and triangular arrays some care is necessary. The 
energy required to store an additional electron on an island is e^/2CE, where Cs is the sum 
of the capacitances to other islands and to ground. As in triangular arrays all islands are 
coupled with z = 6 instead of z = 4 junctions, the required energy is 2/3 times smaller than 
that of an island in a square lattice. Similarly, the freedom of the phase on a particular 
island is determined by the Josephson coupling energy of all junctions connected to the 
island and therefore it seems reasonable to assume that in a triangular array the effective 
Josephson coupling energy is 3/2 times that of a square array. Summarizing, we come to the 
conclusion that the effective Ec/ Ej-Ta.tio for a triangular array is a factor 4/9 lower than 
that of a square array. 



APPENDIX C: PHASE CORRELATOR 



In this Appendix we show how to evaluate the phase-phase correlator gij introduced in 
Section [11 B| (see also Ref. [ p.21|| ). The starting point is: 



gio{r) = (exp[0i(r) - 0o(O))ch 



By making use of the parametrization 

T 

(/)i(r) = (/)jo + 27rmj- + 6'j(r) 

it is possible to verify that all the off-diagonal elements of the correlation function, viz. 
giQ{j) for 2 7^ vanish because of the integrations over 0jo- The reason is that Sch does not 
depend on the phase 0o(''") itself but only on its time derivative. It is therefore sufficient to 
calculate the on-site correlation function at site 

1 / 47r^ \ 

gij) = gooi-r) = — XI exp (2m ^ q^jUj - exp(-27riTnor) gdr) (C2) 

^c'* {nj} j ij 
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where gdT), the correlation function for the case of continuous charges, results from the 
remaining integral over 6{t). It is given by 

g^r) = exp [ - 2e'CoM^ - rT)] . (C3) 

By using the Poisson resummation formula 



{nj ij i ^ ^^^^ {gj ij 

and performing the Fourier transform 

g{uju) = / dTexp{iUi,T)g{T) 







one obtains 



The expressions for the coefficients e, 7, C and A in the coarse-graining approach follow from 
an expansion at small frequencies of Eq. ( P5|) . 



APPENDIX D: DERIVATION OF THE COUPLED COULOMB GAS ACTION 

In this Appendix we briefly discuss the steps leading to Eq.(^). following the calculation 
in Ref. [pT] , p2[ . First a path integral representation is introduced for the island charges. In 
terms of phase trajectories 0j(T) and charges qi{T) = Qj(r)/2e the partition function takes 
the form 

^ = n / dqjo J Dq, UT. J DH^) exp(-5{0, g}), (Dl) 

j i {m,} 

where the phases obey the boundary conditions 0j(O) = + 27mj, while the charge 

paths are periodic, gj(0) = qj{l3) = qjo. The effective action in the mixed representation is 



S{4>, q} = r rfr I 2e' q^{T)C;^ q,{r) + i J] g.(r)0,,(r) - J] cos(0, - 0,) I . (D2) 

° I id * {ij) J 



Summation over winding numbers {rii} fixes the charges to be integer-valued. Starting 
from the partition function ( |D1| ), we first introduce the vortex degrees of freedom. This 
can be done by means of the Villain transformation (see also |^); the time-dependent 
quantum problem requires some additional steps We introduce the lattice with 



spacing e in time direction; this spacing is of order of inverse Josephson frequency: e ^ 
{SEjEqY^'^. In the Villain approximation one replaces 

exp { - eEj cos(0i,^ - 0j,r)]} ^ X! ^^P { — — ~ 27rmi^|^}. 



(D3) 
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Here, we have introduced a two-dimensional vector field iiii^, defined on the dual lattice 
(alternatively, it can be considered as a scalar field defined on bonds). The function 

1 1 



Fix) 



a; < 1, 



2xln(Jo(x)/Ji(x)) 2xln(4/x) 

has been introduced to "correct" the Villain transformation for small Ej. As we see, its 
entire effect is to renormalize (increase) the Josephson coupling Ej EjF{eEj), but it 
does not affect the physics. In the following we will use only the renormalized constant. 
The rhs. of Eq. (03) can be rewritten as 



2eEj 



i3irV(f)ir] 



Now the Gaussian integration over the phases can be easily performed, yielding 

liT Jir i,j,T J i,T 

and the summation is constrained by the continuity equation, 

V ■ Jj^ - = . 

The time derivative stands for a discrete derivative /(r) 
constraint is satisfied by the parameterization |^ 

4^)=n('^)(nV)-ig,, + e('''^)V,A,.. 



(D4) 



^-Af{r + e,)-f{r)]. The 



Here the operator (nV)~^ is the line integral on the lattice (in Fourier space it has the form 
i{kx + ky)~^)^ e^^'^'^ is the antisymmetric tensor, while Aij- is an unconstrained integer-valued 
scalar field. It is important to mention that the continuity equation can be solved in different 
ways which can be mapped onto each other by gauge transformations (see Ref. ||90|| ). 

With the use of the Poisson resummation (which requires introducing a new integer 
scalar field Vir) the partition function can be rewritten as 

Z= exp-S{q,v}. 

[giT,'UiT] 

The effective action for the integer charges and vorticities Vi is 



S[q,v] = 2e'eY^q,^Cr^q,^ " ^E KHnV)-^^ 



47r 



ijT 



TTVi, 



V,n^^\nV)-'q^. 



G 



eE 



The kernel Gij is the lattice Green's function, i.e., the Fourier transform of k~'^. Finally, 
after some algebra one arrives to the effective action of Eq.(^), which we rewrote, for 
simplicity, in the continuous notations. 
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APPENDIX E: EFFECTIVE SINGLE VORTEX ACTION 



It is possible to derive from the coupled Coulomb gas action an effective action for a 
single vortex of vorticity v = ±1. The single vortex effective action includes the effect of 
the interaction with fluctuating charges and the other vortices which are present in the 
system due to quantum fluctuations. The desidered effective action is formally obtained 
by performing the sum in the partition function over all charge and vortex configurations 
excluding the vortex whose dynamics is to be studied. We introduce the vortex trajectory 

Vi^r = vS{ri - r{T)) . 

The single vortex effective action can then be written as 

Seff = - In (6271 Ejv ^ Vi^rGij6{rj - r(r)) + 

+ivY^qi^rQiA^]-r{T)) + iveY,\r-VQij8{rj-r{T))) , (El) 

where the average is to be taken with the full coupled Coulomb gas action. The first 
term describes the static interaction with other vortices, whereas the second describes the 
dynamical interaction with charges. This expression, although exact, cannot be evaluated 
explicitly because of the nonlinearity of the action. To proceed we expand the dynamical 
part of the average in Eg. ([El|) in cumulants up to second order in the vortex velocity r(r). 
For a uniform external current distribution the result is 

Seff = ^ E nr)MaMr) - r(r'), r - r')r\r') + 2mve ^ e„,/>^(r) , 



rr 



r 



Mab = E Vae(r(r) - r,){q^,qkr')VbQ{rk - r(r')) , (E2) 

jk 

where a,b = x,y and eat is the anti-symmetric tensor |pi5 |. 
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APPENDIX F: LIST OF SYMBOLS 



Josephson Energy Ej 

Junction capacitance C 

Ground capacitance Cq 
Capacitance matrix 

Charging energy (junction) Ec = e^/(2C) 

Charging energy (ground) Eq = c^Cqq /2 

Index labels for the islands in the array i,j 

Superconducting order parameter Ae"^ 

Charge on the island Q 

External charge Qx 

Vector potential A^j = // A • dl 

Flux quantum $o — h/2e 

Magnetic frustration per plaquette / 

External current / 

Junction critical current /c 

Quantum of resistance it^g = /i/(4e^) 

Dissipation strength a 

Josephson plasma frequency cUp — yf&EjEc 

BCS transition temperature 

Vortex unbinding transition Tj 

Charge unbinding transition Tch 

Universal conductance at the S-I transition a* 

Vortex mass My 

McCumber parameter /3c 

Josephson Junctionn Array JJA 

Superconductor - Insulator S-I 

Berezinskii - Kosterlitz - Thouless BKT 

Aharonov - Bohm AB 

Aharonov - Casher AC 

Bose - Hubbard BH 

Quantum Phase Model QPM 
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FIGURES 




Fig 1 - R.. F^jrio ajid. H. van do" Zant. 



FIG. 1. Left hand side: a Scanning-Electron-Microscope (SEM) photograph of the Heisenberg 
transistor. The leads in upper, left corner are used to perform a four-terminal measurement 
on the two junctions in series. Offset-charges on the central island are nulled out by the gate 
capacitor situated below the central island. (Picture taken by W.J. Elion.) Right hand side: 
a schematic drawing of the Heisenberg transistor ||T^. The phase and the charge on the island 
are quantum mechanical conjugated variables. By varying the flux through the SQUID ring the 
effective Josephson energy is tuned and, as a result, phase fluctuations on the central island can 
be varied. (Reprinted by permission from Nature 371, 594 (1994) copyright 1994 Macmillan 
Magazines Ltd.) 
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FIG. 2. A Scanning-Electron-Microscope (SEM) photograph of a two-dimensional, square array 
of small tunnel junctions produced by shadow evaporation. The white bar is 1 /xm long. (Picture 
taken from the Ph.D. thesis of L.J. Geerligs, Delft 1990 (unpublished).) 
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FIG. 3. A Josephson array consists of a regular network of superconducting islands weakly 
coupled by tunnel junctions. Each junction is characterized by the Josephson coupling Ej and the 
junction capacitance C; each island by the capacitance to the ground Cq. 
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Fig 4 - R. Fasio ajid. H. van der Zajut 



FIG. 4. The zero-field linear resistance per square measured as a function of the temperature 
for six different square arrays. The curved dashed lines are fits to the vortex-BKT square root cusp 
formula. The dashed horizontal line indicates the zero-temperature universal resistance at the S-I 
transition calculated in Ref. |15[l(| . (From Ref. []5l| .) 
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Fig, 5 - R. Fazio and H. van der Zant 



FIG. 5. Measured phase diagram in zero magnetic field for square (solid squares) and triangular 
(solid triangles) arrays showing the S-I transition at Eq/Ej ~ 1.7. The solid line is a guide to 
the eye connecting the data points and the dotted line on the superconducting side represents the 
result of the calculation of Ref. ||76|. (From Ref. [pl|].) 
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FIG. 6. Sketch of the phase diagram for short-range interaction between charges. The phase 
incoherent state is resistive with an activated behavior of the resistance. At T = the array is a 
insulator. The dependence of the quantum critical point on the capacitance matrix is all contained 
in £^0- 
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FIG. 7. The phase diagram for a quantum JJA in the hmit of long-range (logarithmic) inter- 
action between charges. Similarly to vortices the charges undergo a BKT transition leading to 
insulating behavior at low temperature and small Ej/Eq ratio. 
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FIG. 8. Zero-bias resistance versus magnetic frustration for a square (a) and a triangular (b) 
array. The dip at / = 1/2 is the most pronounced one in both figures. In the triangular array 
the dip at / = 1/4 is more pronounced than the one at / = 1/3. In the square array the opposite 
occurs. (From Ref. []5T|.) 
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Fig 9 - R. Fasio and H. van d^r Zant 



FIG. 9. Phase boundaries between the Mott insulating phase (below each curve) and the su- 
perconducting phase (above). Boundaries for various ratios of the junction capacitance C to the 
self-capacitance Co are shown: C/Cq = O.OOOl(n), O.l(o), 0.2(A), and l.O(O). (From Ref. [ pTof .) 
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fig 10 - R Fazio ajid H. van der Zajit 



FIG. 10. Measured phase diagram for square arrays at / = 1/2, showing the S-I transition at 
Ec/Ej ~ 1.2. In the figure the temperature axis is normahzed to the Josephson coupling, i.e., 
r = T/Ej. The sohd and dashed lines are guides to the eye connecting the data points. (From 
Ref. |1|.) 
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Fig 11 - R Fasio and. H. van der Zant 



FIG. 11. Measured phase diagram for square arrays in a magnetic field. Below the dotted line 
samples become superconducting at low temperatures; above this line samples become insulating. 
At non-commensurate magnetic fields, the S-I transition is not sharp and there is an additional 
(intermediate) metallic region not shown in the figure. (From Ref. |51|.) 
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FIG. 12. The T = phase diagram in the hmit of on-site interaction as a function of the charge 
frustration. At the values of Qx for which two charge states are degenerate, the superconducting 
phase extends to arbitrary small Josephson coupling. 
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FIG. 13. The T = phase diagram calculated with on-site interaction and a small nearest 
neighbor charging term Ei. Around qx = 1/2 the half-integer lobe appears. Inside the lobes, 
each number represents the number of Cooper pairs on a particular island. For example, the 
intermediate lobe, centered around Qx = 1/2 has a checkerboard structure. 
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Fig 14 - R. Fazio ajid H. van der Zajit 



FIG. 14. Mean-field phase diagram in the presence of charge frustration obtained in the 
hard-core limit. As discussed in the text the fictitious field h is related to Qx — 1/2 and the 
coupling J corresponds to Ej. In the figure nn and nnn charging terms are different from zero 
E2 = O.lEi. The central lobe corresponds to the half- filling case and the two small ones on the 
side are the quarter- filling lobes. Finally SSI and SS2 are different types of supersolid. (From 
Ref. H.) 
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Fig 15 - R Fazio and H. van der Zant 



FIG. 15. Top: Mean-field phase diagram for soft-core bosons, as obtained from the analysis 
of the QPM with on-site and nearest-neighbor (Ei/Eq = 1/5) interaction. The symbols are the 
Monte Carlo data (the region between the dotted circles and the crossed circles is the supersolid). 
The checkerboard charge-density wave is denoted by "Sol", the supersolid phase by "SSol", the 
superfluid phase is denoted by "SF" and the Mott-insulating phase by "MI" . Bottom: Supersolid 
region "SSol" at qx= 0.5 as a function of Ui/Uq in the mean-field approximation of. Inset: Oc- 
cupation-number probability |c2p at qx= 0.5 for the two sub-lattices A and B at the particular 
value of Ei/Eq = 0.2. (The notation is slightly different from that used in this review, Uq — > Eq, 
Ui -^Ei,J^ Ej, no ^ (Zx). (From Ref. B.) 
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FIG. 16. Measured low-temperature phase diagram for fifteen square arrays plotted as the 
Ej/ Ec-'catio vs. Rq/Rn- Samples ^ 1 — 8 show a decreasing resistance as temperature is lowered 
(superconducting arrays). Samples jj 9 — 15 show an increasing resistance as temperature is lowered 
(insulating arrays). This classification roughly agrees with the theoretical prediction |91| which 
states that insulating arrays are to be found in the region bounded by the dotted line. Only the 
insulating samples [J 9 — 11 fall slightly outside this area. The diagonal line represents Aq = 2£'c 
and the dashed line corresponds to a Stewart McCumber parameter of (5c = 1- (From Ref. |161|.) 
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Fig 17 - R Fasio ajid H. van der Zant 



FIG. 17. The phase diagram at T = in the local damping model. Along the solid line the 
conductivity is universal, whereas it is a function of the dissipation strength along the dotted line. 
(From Ref. |7|.) 



101 



-30 



-2C 



-1& 10 



30 



30 



Fig 18 - R Fasio and. H, van der Zaiit 



FIG. 18. Measured array I — V characteristics for eight values of the back-gate voltage corre- 
sponding to different ground plane resistances Rg. The I — V .s change from superconducting-like 
to insulator-like function oiRg. (From Ref. [|153| .) 
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Fig 19 - R Fasio and. H, van der Zaiit 



FIG. 19. Measured zero-bias array resistance Rq as a function of the ground plane resistance 



Rg. The inset shows the temperature dependence of the resistance. (From Ref. |153|. 
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Fig 20 - R Fasio and. H, van der Zaiit 



FIG. 20. Calculated phase diagram of an array coupled capacitively to a 2DEG. The insets 
show the resistance as a function of the temperature in the different regions of the phase diagram. 
(From Ref. [|l7|.) 



104 



Fig 21 - R Fasio and. H, van der Zaiit 



FIG. 21. A Scanning-Electron- Microscope (SEM) image of a SQUID chain. (From Ref. [I97|). 
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Fig 22 - R Fasio and. H, van der Zaiit 



FIG. 22. Calculated phase diagram of a Josephson chain shown as a function of the Josephson 
coupling and the dissipation strength. The dissipative part of the action leads to two new phases. 
For a > 1/2 the dipoles are bound in a gas of quadrupoles and moreover for strong dissipation 
there is an additional phase transition which separate the quadrupole phase from a phase in which 
the system shows local order. (From Ref. [204|). 
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Fig 23 - R. Fasio ajid H. van der Zant 



FIG. 23. Resistance vs. temperature for eight values of the magnetic field in the range — 64 G. 
The two sets of curves correspond to two different one-dimensional arrays with N = 255 (solid lines) 
and = 63 (dashed lines) junctions. The longer array shows a sharper S-I transition. At the 
point J* (J = y^Ej/Ec) the resistance is length independent. (From Ref. |197[| .) 
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Fig 24 - R Fasio and H. van der Zant 



FIG. 24. The magnetic field dependence of the threshold voltage for one-dimensional arrays of 
different length. The inset shows the corresponding I — V curves. (From Ref. |197|j.) 
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Fig 25 - R Fasio aiLd H. van der Zant 



FIG. 25. Temperature dependence of the resistance for a square (a) and triangular (b) 2D 
Josephson array measured at different applied magnetic fields. The field tuned S-I transition 
occurs at that frustration where the temperature dependence of the resistance changes sign. In 
both cases this change occurs at / = 0.10 — 0.15. (From Ref. |5l||.) 
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Fig 26 - R. Faado and. H. van der Zant 



FIG. 26. I — V characteristics measured at low temperature (10 mK) for three values of the 
applied field. The square 2D array has an {Ec/Ej ratio of 1.25. The crossover from the su- 
perconducting to the insulating behavior is related to the field tuned S-I transition. (From Ref. 
51.) 
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Fig 27 - R Fasio and. H. van der Zant 



FIG. 27. Field dependence of the resistance of a triangular 2D array measured at different 
temperatures, T = 50 mK (solid line), T = 120 mK (dotted line), and T = 160 mK (dashed line). 
The field-tuned transition is observed around different fractional values of the frustration indicated 
by the open circles. (From Ref. |]5l|| .) 
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Fig 23 - R. Fasio ajid H. van der Zant 



FIG. 28. Resistance as a function of the scaling parameter \f — f^\T~^^^^'^^ for two different 
arrays (data taken in the range < / < 0.2). The data collapse onto a single curve: the upper 
part for the insulating transition and the lower part for the superconducting transition. The inset 
in the figure on the left shows the scaling for the array close to full frustration (data taken in the 
range 0.5 < / < 0.6). The inset in the figure on the right shows a log- log plot of the critical 
frustrations fc as a function of the BKT transition temperature for the four measured samples. 
The line through the data yields a critical exponent z = 1.05. (From Ref. fs^ .) 
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FIG. 29. Phase configuration of a vortex configuration. The arrows indicate the phase of each 
island with respect to a given reference direction. 
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Fig 30 - R Fasio and. H. van der Zant 



FIG. 30. A current-voltage characteristic of an underdamped 2D square, aluminum array mea- 
sured at low temperature (10 mK) in a magnetic field of 0.1 <I>o applied per cell (/ = 0.1). The 
arrow at the left indicates the expected depinning current of Q.lNIc with the number of junc- 
tions perpendicular to the direction of the current flow. For small voltages hysteresis is seen. The 
flux-flow region is found above the depinning current but below the current at which row switching 
sets in (arrow at the left). (From Ref. 236 |.) 
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Fig 31 - R Fasio and H. van der Zant 



FIG. 31. Longitudinal resistance Roxx (a) and the Hall resistance Roxy (b), and the Hall angle 
Q (c) as a function of frustration. Rqxx, and Roxy are shown for various temperatures ranging from 
T = 20 (top), 75, 100, 125, 150, 175 mK. Rqxx is symmetric around / = and / = ±1/2 whereas 
Roxy changes sign upon passing through these frustrations. (From Ref. [161| .) 
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Fig 33 - R. Fasio and H. van der Zant 



FIG. 32. Sample lay-out used to measure ballistic vortices in 2D Josephson arrays. (From 
Ref. @.) 
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Fig 33 - R Fasio and H. van der Zant 



FIG. 33. Voltage vs. Temperature measured across the channel (a), between probes V3 and 
VIO (b), between probes V5 and V8 (c), and between probes V3 and V4 (d). In the inset voltages 
are plotted in units of the voltage across the channel. The dashed line corresponds to the voltage 
across V3-V10 when the current direction is reversed. At low temperatures, the voltage across 
the two probes opposite from the channel is almost equal to the voltage across the channel: all 
vortices that go through the channel leave the array between V7 and V8. With reversed current 
direction vortices are accelerated in the opposite direction and no ballistic motion is observed. 
(From Ref. Q.) 
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FIG. 34. Measured zero-bias resistance per junction versus the inverse normalized temperature 
measured for two different square arrays. At low temperatures the resistance of the sample with the 
smaller Ej/Ec ratio (b) is temperature independent indicative for quantum tunneling of vortices. 
(From Ref. |3|.) 
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Fig 3S - R Fasio and. H, van der Zaiit 



FIG. 35. Schematic drawing of the hexagon-shaped Josephson array to measure vortex inter- 
ference. (From Ref. [p66|.) 
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Fig 36 - R Fasio aiLd H. van der Zant 



FIG. 36. A sketch of the sample layout of a quasi one-dimensional Josephson array. The current 
is injected in the middle while the voltage probes are situated at the end of the busbars. (From 
Ref. |6|.) 
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Fig 37 - R. Fasio and H. van der Zant 



FIG. 37. Schematic drawing of the energy bands for a vortex moving in a quasi-lD Josephson 
array. Dots: numerical calculated energy bands starting from Schrodingers equation with a cosine 
potential. The dashed line shows the first band of a cosinusoidal dispersion relation with the same 
band width. (From Ref. |23(]f| .) 



121 



.V tMV> 



iff 38 ^ R f 'asio Bud H, yaji dgr 2ajit. 



FIG. 38. A nose-shaped I-V characteristic indicative for Bloch oscillating vortices. Data (cir- 
cles) has been obtained in a quasi-lD array with size 7 by 1000 and has been measured at 10 mK 
for a ID vortex density (n) of 0.04. The solid line is the analytical result discussed in the text. 
(From the Ph.D. thesis of A. van Oudenaarden, Delft, 1998, unpublished.) 
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FIG. 39. Zero-bias resistance Rq as a function of the one-dimensional vortex density (n = Wf) 
for two samples with different widths W. The sample length is 1000 cells and data has been 
obtained at 30 mK. The bottom curve (W = 7) is mirrored with respect to the x-axis for clarity. 
(From Ref. [p65| .) 
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Fig 4fl - R Fasio and. H, van der Zaiit 



FIG. 40. Activation energy as a function of n — no {uq is the commensurate density) in the 
case that the boundary pinning Eb dominates over the sohton formation energy Es- On the 
incommensurate side of the transition, n > nc, sohtons form spontaneously and the physics is 
determined by boundary pinning and the elastic energy. (From Ref. |286(| .) 
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Fig 41 - R Fasio and. H, van der Zaiit 



FIG. 41. Activation energy as a function of n — uq when the sohton formation energy 
£Jsdominates over the boundary pinning Ef, (opposite hmit as considered in Fig. |40| ). (From 
Ref. M.) 
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FIG. 42. Measured activation energy of an array of 1000 x 7 cells with Ej = 0.9 K and Eq = 0.7 
K. The dashed line is a fit to the data yielding the width of the Mott region. The inset shows En 
inside the Mott phase. (From Ref. [p86| .) 



126 



T(mK) 
flOO 200 LOO 



1 ! 






- \ 


c i = 




D 


^ 6=1 
=■ 5=2 




|. B F. 




H - 


■ . 1 , L 




1 ■ 


2 + 


{, 





E/kjT 



Ftg 43 - R Fasio and H. van der Zant 



FIG. 43. Arrenius plot of the linear resistance for the ordered (squares) and for the disordered 
(triangles, circles) arrays. At low temperatures (right hand side of the figure), the resistance of 
the disordered arrays (triangles and circles) has dropped below measuring accuracy (dashed line) 
indicating vortex localization. (From Ref. |264| .) 
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Fig 44 - R Fario and. H. van der Zant 



FIG. 44. The experimental setup to detect interference effects of vortices: The Corbino disk 
(FVom Ref. Q.) 
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FIG. 45. a) A charge qubit. b) Improved qubit design as proposed by the Karlsruhe group. 
The island is coupled to the circuit via two Josephson junctions with parameters Cj and Ej. This 
dc-SQUID is tuned by the external flux which is controlled by the current through the inductor 
loop (dashed line) . The setup allows switching the effective Josephson coupling to zero. (Reprinted 
by permission from Nature 398, 305 (1999 copyright 1999 Macmillan Magazines Ltd.) 
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Fig 45 - R Fasio and. H, van der Zaiit 



FIG. 46. The three-junction flux qubit. Josephson junctions 1 and 2 have the same Josephson 
energies Ej and capacitance C; Josephson junction 3 has a Josephson energy and capacitance that 
are a times larger. The islands are coupled by gate capacitors Cg = 7C to gate voltages Va and 
Vb- The arrows define the direction of the currents. (From Ref. ||305|| .) 
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